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ABSTRACT 

We examine a broad range of mass flow diagnostics in a large sample of Berbig Ae/Be stars 
(BAEBES) using high resolution optical spectra. The B/3 and Be I 5876 A lines show the high¬ 
est incidence of P-Cygni (30%) and inverse P-Cygni (14%) morphologies, respectively. The Fe B 
4924 A line also shows a large incidence of P-Cygni profiles (11%). We find support for many of 
the conclusions reached in a study based on the analysis of the Be I A10830 line in a large sample 
of BAEBES. Namely, BAEBES exhibit smaller fractions of both blue-shifted absorption (i.e. mass 
outflow) and red-shifted absorption (i.e. mass infall or accretion) than their lower mass cousins, the 
classical T Tauri stars (CTTSs). In particular, the optical data supports the conclusion that BAEBES 
displaying red-shifted absorption, in general, show maximum red-shifted absorption velocities that 
are smaller fractions of their stellar escape velocities than is found for CTTSs. This suggests that 
BAEBE accretion flows are originating deeper in the gravitational potentials of their stars than in 
CTTS systems. In addition, we find a lack of inner disk wind signatures in the blue-shifted ab¬ 
sorption objects; only stellar wind signatures are clearly observed. These findings, along with the 
lack of detected magnetic fields around BAEBES, support the idea that large magnetospheres are 
not prevalent around BAEBES and that accretion flows are instead mediated by significantly smaller 
magnetospheres with relatively smaller truncation radii (e.g. 1-2 i?*). Red-shifted absorption is much 
more common around Berbig Ae stars than Be stars, suggesting that Berbig Be stars may accrete via 
a boundary layer rather than along magnetic field lines. 

Subject headings: accretion-stars:pre-main sequence-stars:variables:T Tauri, Berbig Ae/Be- 
stars:winds,outflows-methods:statistical 


1. INTRODUCTION 

The interplay between mass loss and accretion in pre- 
main sequence systems is important in determining both 
the final stellar mass, the mass in the disk available to 
form planets, and the local physical conditions (e.g. in¬ 
cident UV flux, temperature) affecting planet formation 
and evolution. While many details of the physical mech¬ 
anisms involved in mass accretion and mass outflow in 
the immediate circumstellar environments (t< 0.1 AU) 
for low-mass Class B pre-main sequence stars, or clas¬ 
sical T Tauri stars (CTTSs), are now understood (e.g. 
Muzerolle et al. 2001; Romanova et al. 2009), the pic¬ 
ture is less clear concerning these processes around the 
intermediate mass Berbig Ae/Be stars (BAEBES). 

Optical spectroscopic investigations into the nature of 
the circumstellar environments around BAEBES have 
been numerous. A small number of these studies have 
been performed on medium or large samples of BAEBES. 
Finkenzeller & Mundt (1984) presented medium resolu¬ 
tion (Av~16 km s“^) Ba and Na I D lines for a sam¬ 
ple of 43 BAEBES that focused mainly on the P-Cygni 
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profiles of a small number (8) of objects. The investi¬ 
gations of Bamann & Persson (1992), Bohm & Catala 
( 1994 ), and Bohm & Catala (1995) all examined sam¬ 
ples of ^30 BAEBES and a small or medium number of 
spectral line diagnostics. Bamann & Persson (1992) and 
Bohm & Catala (1994) compare their samples to known 
TTS characteristics but these studies focus mainly on 
pure emission line morphologies and not on red and blue- 
shifted absorption. Vieira et al. (2003), in a spectroscopic 
BAEBE study covering the largest number of objects to 
date, classified the Ba profiles of 131 stars using low and 
medium resolution spectra. Vieira et al. (2003) also ex¬ 
amined forbidden line emission in their sample. Bowever, 
no detailed interpretation of the Ba profile morphologies 
or forbidden line strengths are performed, yielding little 
information about the kinematics of the inner circum¬ 
stellar material. Finally, the BAEBE variability study 
of Mendigutia et al. (2011) examined Ba, [O I]6300, Be 
I 5876, and the Na I D lines in a medium sized sample us¬ 
ing medium resolution spectra. Bigh spectral resolution 
observations of a wide range of accretion and outflow di¬ 
agnostics in a large sample of BAEBES is clearly needed. 

Besides the studies listed above, most optical spectro¬ 
scopic BAEBE studies have focused on single objects or 
small sample sizes and a small number of spectral line 
diagnostics (e.g. Catala et al. 1993; Kraus et al. 2008; 
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Grady et al. 2010). While this is useful for providing de¬ 
tails about specific objects, studies of this sort provide lit¬ 
tle information concerning the physical characteristics of 
mass flows around HAEBES as a group, which are impor¬ 
tant for illuminating the general differences between the 
early evolution of low and intermediate mass pre-main 
sequence stars. Comparing large sample of HAEBES and 
CTTSs can also provide insight into the important phys¬ 
ical mechanisms responsible for generating mass flows in 
young stars. 

While both HAEBES and CTTSs host gaseous accre¬ 
tion disks, the lack of strong detected magnetic fields 
on HAEBES (e.g. Wade et al. 2007; Alecian et al. 2013) 
and their large usini values compared to CTTSs Johns- 
Krull (2007) provide hints that accretion and outflow 
processes do not proceed identically in both groups of 
objects. Magnetospheric accretion, widely accepted to 
operate in most CTTS systems, funnels material from a 
truncated inner disk onto the surface of the star (Konigl 
1991). Accretion along magnetic field lines is not possible 
from outside of a narrow region around the corotation ra¬ 
dius in the disk (Shu et al. 1994). Muzerolle et al. (2004) 
pointed out that in order for magnetospheric accretion 
to operate around HAEBES, the large nsini values of 
HAEBES require the corotation radii, and thus the mag¬ 
netic truncation radii, to be closer to the stellar surface 
than in CTTSs. This has implications for the ability of 
HAEBES to launch accretion-generated outflows in the 
same way as CTTSs (see §of Cauley & Johns-Krull 2014, 
for an expanded discussion). Spectral line profile mor¬ 
phologies can provide insight into how these mass flows 
are generated. 

In this paper we present a high resolution optical study 
of a large sample of HAEBES using a large number of 
spectral line diagnostics with the goal of investigating dif¬ 
ferences in the accretion and outflow mechanisms oper¬ 
ating around HAEBES compared to CTTSs. This work 
is the second part of a multi-wavelength study of ac¬ 
cretion and outflows in HAEBES. The first part is pre¬ 
sented in Cauley & Johns-Krull (2014, hereafter Paper 
I) , and focused on the study of He I A10830 in a sample 
of 56 HAEBES. Our comprehensive coverage of optical 


mass flow indicators enables us to derive reliable esti¬ 
mates and confidence intervals of the incidence of accre¬ 
tion and outflow signatures in HAEBES. We describe our 
observations and data reduction in §2. In §3 we discuss 
the calculation of stellar radial velocities and rotational 
velocities. Section 4 covers the classification of the line 
profiles. We present the individual line and overall mass 
flow statistics in §5, as well as a comparison to the He I 
A10830 statistics from Paper I. A comparison of the en¬ 
tire sample statistics to those of CTTSs is also be given 
in §5. Section 6 investigates the physical cause of the 
differences in morphology statistics. Our conclusions are 
summarized in §7. 

2. OBSERVATIONS AND DATA REDUCTION 

Our sample consists of 87 HAEBES covering a wide 
range of spectral types and was chosen from the cata¬ 
logues of Vieira et al. (2003); The, de Winter, & Perez 
(1994); and Finkenzeller & Mundt (1984). Classihcation 
as a HAEBE in a previous survey was the only selec¬ 
tion criteria. It is important to note that our sample 
likely contains objects at various phases in their pre- 
main sequence evolution. However, the small number of 
HAEBES currently identified in the literature weakens 
the signihcance of a more specific comparison between 
distinct evolutionary groups. Thus we include all ob¬ 
jects in our central analysis regardless of spectral type 
and estimated age, though presumably all are pre-main 
sequence in nature. 

2 .1. Observations 

Our observations were performed using two cross- 
dispersed echelle spectrographs on two telescopes: the 
Sandiford Echelle Spectrometer (SES McCarthy et al. 
1993) on the McDonald Observatory 2.1 m Otto Struve 
telescope and the Tull 2dCoude echelle spectrometer 
(TS2 Tull et al. 1995) on the 2.7 m Harlan J. Smith 
telescope. The specifics of each observation are given 
in Table 1. Some objects were observed multiple times 
during different observing runs. The use of observations 
from different epochs in the analysis is described in sec¬ 
tion 5. 


Table 1 

Log of echelle observations 


Object ID 
(1) 

Instrument 

(2) 

Telescope 

(3) 

UT Date 
(4) 

Integration time (s) 

(5) 

S/N @ 6400 A“ 
(6) 

AB Aur 

TS2 

2.7 m 

18-Sep-2011 

120 

65 

AE Lep 

SES 

2.1 m 

Ol-Nov-2012 

2400 

35 

BD-l-61 154 

TS2 

2.7 m 

20-Sep-2011 

2500 

45 

» 

SES 

2.1 m 

30-0ct-2012 

1200 

30 

55 

SES 

2.1 m 

31-Oct-2012 

2000 

40 

BE Ori 

TS2 

2.7 m 

21-Jan-2013 

3000 

65 

BH Cep 

TS2 

2.7 m 

20-Sep-2011 

2500 

45 

!) 

SES 

2.1 m 

Ol-Nov-2012 

3600 

30 

5 ? 

SES 

2.1 m 

O 2 -N 0 V- 2 OI 2 

3600 

30 

CQ Tau 

TS2 

2.7 m 

19-Sep-2011 

2100 

65 

DW CMa 

TS2 

2.7 m 

20—Jan—2013 

3600 

20 

CSC 04794-00827 

TS2 

2.7 m 

21-Jan-2013 

3600 

50 

HD 141569 

TS2 

2.7 m 

20-Sep-2011 

180 

90 

!) 

TS2 

2.7 m 

20—Jan—2013 

75 

60 

HD 142666 

TS2 

2.7 m 

20—Jan—2013 

600 

65 

HD 163296 

TS2 

2.7 m 

18-Sep-2011 

300 

55 

HD 169142 

TS2 

2.7 m 

19-Sep-2011 

200 

50 

HD 190073 

TS2 

2.7 m 

18-Sep-2011 

180 

55 

5 ) 

TS2 

2.7 m 

26-Sep-2012 

350 

60 

HD 203024 

TS2 

2.7 m 

20-Sep-2011 

390 

60 
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Table 1 — Continued 


Object ID 
(1) 

Instrument 

(2) 

Telescope 

(3) 

UT Date 
(4) 

Integration time (s) 

(5) 

S/N @ 6400 A“ 
(6) 


SES 

2.1 m 

30-0ct-2012 

600 

30 


SES 

2.1 m 

31-Oct-2012 

1500 

50 

HD 244314 

TS2 

2.7 m 

20-Sep-2011 

2200 

65 


TS2 

2.7 m 

20—Jan—2013 

3000 

80 

HD 244604 

TS2 

2.7 m 

20-Sep-2011 

1300 

75 

HD 245185 

TS2 

2.7 m 

19-Sep-2011 

2300 

65 

HD 249879 

TS2 

2.7 m 

26-Sep-2012 

3200 

60 

” 

SES 

2.1 m 

Ol-Nov-2012 

3600 

35 


SES 

2.1 m 

O 2 -N 0 V- 2 OI 2 

3600 

45 

HD 250550 

SES 

2.1 m 

Ol-Nov-2012 

2000 

40 

” 

SES 

2.1 m 

O 2 -N 0 V- 2 OI 2 

3000 

100 

HD 287823 

TS2 

2.7 m 

22-Sep-2011 

3000 

55 

HD 290409 

TS2 

2.7 m 

22-Sep-2011 

3000 

70 

HD 290500 

TS2 

2.7 m 

23-Sep-2011 

3600 

60 

HD 290764 

TS2 

2.7 m 

23-Sep-2011 

2400 

70 

HD 290770 

TS2 

2.7 m 

21-Sep-2011 

1200 

90 

HD 34282 

TS2 

2.7 m 

21-Sep-2011 

1500 

70 

HD 35187 

TS2 

2.7 m 

18-Sep-2011 

240 

70 


SES 

2.1 m 

Ol-Nov-2012 

1200 

35 

?? 

SES 

2.1 m 

O 2 -N 0 V- 2 OI 2 

2400 

40 

HD 35929 

TS2 

2.7 m 

18-Sep-2011 

400 

50 

HD 36408 

TS2 

2.7 m 

18-Sep-2011 

120 

40 

HD 37357 

TS2 

2.7 m 

18-Sep-2011 

500 

55 

HD 37411 

SES 

2.1 m 

30-0ct-2012 

2400 

30 

HD 38120 

SES 

2.1 m 

30-0ct-2012 

600 

30 


SES 

2.1 m 

31-Oct-2012 

2000 

55 

HD 50083 

TS2 

2.7 m 

18-Sep-2011 

180 

50 

HK Ori 

TS2 

2.7 m 

21-Sep-2011 

3600 

30 

IL Cep 

TS2 

2.7 m 

20-Sep-2011 

1300 

80 


TS2 

2.7 m 

22-Sep-2011 

1200 

75 

?? 

TS2 

2.7 m 

23-Sep-2011 

1300 

80 

IP Per 

TS2 

2.7 m 

19-Sep-2011 

2400 

70 

IRAS 05044-0325 

TS2 

2.7 m 

20—Jan—2013 

3200 

35 

IRAS 06071+2925 

TS2 

2.7 m 

20—Jan—2013 

3600 

30 

IRAS 07061-0414 

TS2 

2.7 m 

21-Jan-2013 

3600 

60 

IRAS 17481-1415 

TS2 

2.7 m 

20-Sep-2011 

3600 

30 

IRAS 18306-0500 

TS2 

2.7 m 

20-Sep-2011 

3600 

30 

IRAS 18454+0250 

TS2 

2.7 m 

19-Sep-2011 

3000 

20 

IRAS 19343+2926 

TS2 

2.7 m 

19-Sep-2011 

3000 

20 


TS2 

2.7 m 

26-Sep-2012 

3000 

15 

V 

SES 

2.1 m 

Ol-Nov-2012 

3600 

10 

?? 

SES 

2.1 m 

O 2 -N 0 V- 2 OI 2 

3600 

10 

LkHo; 134 

TS2 

2.7 m 

20-Sep-2011 

2700 

60 


SES 

2.1 m 

Ol-Nov-2012 

3600 

35 


SES 

2.1 m 

O 2 -N 0 V- 2 OI 2 

3600 

40 

LkHa 208 

SES 

2.1 m 

Ol-nov-2012 

3600 

35 


SES 

2.1 m 

O 2 -N 0 V- 2 OI 2 

3600 

35 


TS2 

2.7 m 

22—Jan—2013 

3000 

65 

LkHa 233 

TS2 

2.7 m 

22-Sep-2011 

3600 

20 


TS2 

2.7 m 

26-Sep-2012 

3600 

20 

LkHa 257 

TS2 

2.7 m 

21-Sep-2011 

3800 

35 

LkHa 324 

TS2 

2.7 m 

20-Sep-2011 

3600 

40 

LkHa 339 

TS2 

2.7 m 

22—Jan—2013 

1960 

20 

MWC 1080 

TS2 

2.7 m 

21-Sep-2011 

2500 

20 


TS2 

2.7 m 

22-Sep-2011 

2500 

20 


TS2 

2.7 m 

23-Sep-2011 

2500 

20 


SES 

2.1 m 

30-0ct-2012 

2000 

20 


SES 

2.1 m 

31-Oct-2012 

3000 

10 

MWC 120 

TS2 

2.7 m 

18-Sep-2011 

400 

80 

MWC 137 

TS2 

2.7 m 

20—Jan—2013 

1800 

15 

MWC 147 

TS2 

2.7 m 

26-Sep-2012 

600 

20 

MWC 300 

TS2 

2.7 m 

18-Sep-2011 

2700 

10 

MWC 361 

TS2 

2.7 m 

18-Sep-2011 

180 

70 

MWC 480 

TS2 

2.7 m 

18-Sep-2011 

180 

50 

MWC 610 

TS2 

2.7 m 

18-Sep-2011 

360 

50 


TS2 

2.7 m 

26-Sep-2012 

900 

60 

MWC 614 

TS2 

2.7 m 

18-Sep-2011 

120 

75 

MWC 758 

TS2 

2.7 m 

18-Sep-2011 

400 

55 

MWC 778 

SES 

2.1 m 

30-0ct-2012 

3600 

10 


SES 

2.1 m 

31-Oct-2012 

4000 

10 

MWC 863 

TS2 

2.7 m 

23-Sep-2011 

500 

50 

NZ Ser 

TS2 

2.7 m 

19-Sep-2011 

3000 

40 

R Mon 

TS2 

2.7 m 

26-Sep-2012 

2000 

20 

RR Tau 

SES 

2.1 m 

30-0ct-2012 

2500 

35 


SES 

2.1 m 

31-Oct-2012 

3600 

35 
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Table 1 — Continued 


Object ID 
(1) 

Instrument 

(2) 

Telescope 

(3) 

UT Date 
(4) 

Integration time (s) 

(5) 

S/N @ 6400 A“ 
(6) 


TS2 

2.7 m 

20—Jan—2013 

3600 

55 

SV Cep 

TS2 

2.7 m 

20-Sep-2011 

2500 

65 

T Ori 

TS2 

2.7 m 

23-Sep-2011 

3000 

70 


SES 

2.1 m 

30-0ct-2012 

1200 

20 


SES 

2.2 m 

31-Oct-2012 

3600 

50 

UX Ori 

TS2 

2.7 m 

22-Sep-2011 

3600 

30 

UY Ori 

SES 

2.1 m 

30-0ct-2012 

3600 

20 

VI185 Tau 

TS2 

2.7 m 

19-Sep-2011 

2400 

70 


SES 

2.1 m 

30-0ct-2012 

2000 

30 


SES 

2.1 m 

31-Oct-2012 

3000 

40 

V1578 Cyg 

TS2 

2.7 m 

20-Sep-2011 

2500 

90 


TS2 

2.7 m 

22-Sep-2011 

2400 

90 


TS2 

2.7 m 

23-Sep-2011 

2400 

80 

V1685 Cyg 

TS2 

2.7 m 

18-Sep-2011 

2500 

25 


TS2 

2.7 m 

26-Sep-2012 

3000 

30 


SES 

2.1 m 

Ol-Nov-2012 

2000 

25 


SES 

2.1 m 

O 2 -N 0 V- 2 OII 

3000 

30 

V1686 Cyg 

TS2 

2.7 m 

19-Sep-2011 

3600 

20 

V1787 Ori 

TS2 

2.7 m 

22-Jan-2013 

3600 

35 

V1818 Ori 

TS2 

2.7 m 

22-Jan-2013 

3600 

40 

V346 Ori 

TS2 

2.7 m 

21-Sep-2011 

2400 

75 

V351 Ori 

TS2 

2.7 m 

18-Sep-2011 

650 

50 

V361 Cep 

TS2 

2.7 m 

18-Sep-2011 

2200 

50 

V373 Cep 

TS2 

2.7 m 

21-Sep-2011 

3500 

40 

V374 Cep 

TS2 

2.7 m 

21-Sep-2011 

2800 

30 


TS2 

2.7 m 

26-Sep-2012 

3000 

50 


SES 

2.1 m 

30-0ct-2012 

2000 

30 


SES 

2.1 m 

31-Oct-2012 

2000 

60 

V380 Ori 

SES 

2.1 m 

30-0ct-2012 

2400 

15 


SES 

2.1 m 

31-Oct-2012 

3600 

15 

V586 Ori 

TS2 

2.7 m 

20-Sep-2011 

1500 

75 

V590 Mon 

TS2 

2.7 m 

21-Jan-2013 

3600 

40 

V599 Ori 

TS2 

2.7 m 

22—Jan—2013 

3200 

30 

V718 Sco 

TS2 

2.7 m 

23-Sep-2011 

600 

50 

V791 Mon 

SES 

2.1 m 

Ol-Nov-2012 

2400 

30 

” 

SES 

2.1 m 

O 2 -N 0 V- 2 OI 2 

3600 

45 

VV Ser 

TS2 

2.7 m 

18-Sep-2011 

3000 

45 


TS2 

2.7 m 

21-Sep-2011 

3000 

40 


TS2 

2.7 m 

22-Sep-2011 

3000 

40 

WW Vul 

TS2 

2.7 m 

18-Sep-2011 

1600 

40 


TS2 

2.7 m 

26-Sep-2012 

2000 

70 


SES 

2.1 m 

30-0ct-2012 

600 

25 


SES 

2.1 m 

31-Oct-2012 

1500 

45 

XY Per 

TS2 

2.7 m 

19-Sep-2011 

2000 

65 

Z CMa 

TS2 

2.7 m 

21-Jan-2013 

1800 

20 


^ Continuum S/N for SES spectra taken on 31—Oct-2012 and 02- 

Nov-2012 is estimated at 4800 A. 

The TS2 data were obtained during three separate runs 
in 2011 September (7 nights), 2012 September (3 nights), 
and 2013 January (3 nights). For the TS2 observations, 
we employed the 1.2” slit with the E2 grating to achieve 
a resolving power of i?^60,000, or a velocity resolution 
of ~5 km s“^. The grating configuration provided wave¬ 
length coverage from approximately 3700-10000 A. Ex¬ 
posure times were limited to 1-hour in order to maximize 
the number of objects observed with at least moderate 
signal-to-noise. A S/N of ^50 per pixel is typical for the 
TS2 data, although this varies from S/N~10 up to ~90 
depending on a combination of the object brightness and 
weather conditions during the exposure. ThAr compari¬ 
son lamp spectra were taken multiple times each night to 
be used for wavelength calibrations. Flat field exposures 
were collected at the beginning of each night. Telluric 
and spectroscopic standards were observed periodically 
throughout the course of each observing run. 

The SES data were collected during a single 4-night 


observing run in 2012 October. We used two differ¬ 
ent grating settings, each one on two separate nights, 
to record both ’’red” (5430-6620 A) and ’’blue” (4375- 
4870 A) wavelengths. This was necessary due to the 
SES’s smaller wavelength coverage (~500-1000 A) per 
exposure than that of TS2. Even using two wavelength 
settings, we were not able to observe all of the spectral 
diagnostics that are available in the TS2 data. Thus 
objects that are only observed with the SES are lack¬ 
ing data for some diagnostics. The resolving power for 
each grating setting is i?~60,000. ThAr lamp spectra 
were taken throughout the night and flat field exposures 
were obtained at the beginning of each night. Telluric 
and spectroscopic standards were observed, although the 
spectroscopic standards obtained using the TS2 are used 
interchangeably due to the almost identical velocity res¬ 
olutions of the two instruments. 

2 .2. Data reduction 

All observations were reduced using custom IDL rou¬ 
tines optimized for extracting echelle spectra. Each ob¬ 
ject exposure is bias subtracted and divided by a normal- 
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ized median flat to remove pixel to pixel variations. The 
spectral orders are located and then fit using a order 
polynomial. Each order is optimally extracted which also 
removes most hot pixels and cosmic rays. Wavelength 
calibrations are performed on the ThAr exposures using 
a two dimensional polynomial fit to all spectral orders 
simultaneously. Hundreds of individual lamp lines in the 
exposures are matched to known atlas wavelengths. This 
procedure results in dispersion solutions that are accu¬ 
rate to ~0.10 pixels, or 0.005 A, across the entire spec¬ 
tral range. Telluric absorption was subtracted out of the 
spectral region near 6300 A in order to make identifi¬ 
cation of the O I A6300 line unambiguous. Due to the 
prominence of the Ha and Na I D lines, telluric absorp¬ 
tion in these regions does not affect the classification of 
the line morphologies since these absorption are small 
perturbations of the true line profile shape. Equivalent 
width calculations of Ha lines include the telluric fea¬ 
tures which results in a maximum additional uncertainty 
of ~0.40 A. 

Some exposures contain strong telluric O I emission 
near 5577, 6300, and 6363 A. We have modified our 
echelle reduction program to locate the night sky emis¬ 
sion above and below the spectral trace of the star and 
subtract it from the actual object spectrum. However, 
this routine was not always successful, especially for ob¬ 
jects with profiles that were not centered on the slit or 
for irregularly shaped profiles (e.g. binaries). Thus some 
profiles contain residual night sky emission (e.g. HD 
244314). This extra emission can be identified by its 
narrow (~10 km s“^) width and is easily distinguished 
from the stellar line profile. 

Approximate continuum normalization in each order 
was performed by dividing each order by the blaze func¬ 
tion of the same order, derived from the reduced spec¬ 
trum of the median flat for that particular grating set¬ 
ting. A linear fit is then performed to remove any resid¬ 
ual slope. Some objects display strong photospheric 
Balmer features which have wings that often extend to 
the ends of the spectral order. The continuum normal¬ 
ization in these cases is likely not exactly the true stellar 
continuum. However, identification of non-photospheric 
features in these lines is not problematic with the use of 
a comparison standard since the chosen continuum level 
does not affect the line shape. The Fe H 4352 A and 
N H 6583 A features often lie on the wings of Hy and 
Ha, respectively. These profiles are normalized to the 
continuum in the Balmer lines. 

3. PHYSICAL PARAMETERS 

In order to correctly separate circumstellar from pho¬ 
tospheric features, it is important to have estimates of 
both the stellar rotational velocity (wsinj) and the radial 
velocity (RV) of the object relative to the barycenter. 
For many of the objects in our sample these values have 
been estimated in the literature. For objects that do 
not have previous estimates, we fit rotationally broad¬ 
ened and velocity-shifted synthetic spectra to a single 
spectral order to obtain RV and vsini estimates. The 
order selected is chosen based on the strongest available 
photospheric lines. This method works well for A-type 
and late B-type objects but fails for the earliest spectral 
types. Typically the strong He I A4471 and Mg H A4481 


lines are used to measure usini for rapidly rotating B- 
type stars (Halbedel 1996); for more slowly rotating B- 
type stars, metal lines near 6400 A and 4500 A are used 
(Fekel 2003). Many of our early type objects, however, 
show no absorption features at any of these wavelengths. 
For a typical S/N spectrum in our sample, this method 
can also fail for stars with vsini > 200 km s“^ because 
the features become too weak to reliably measure. Lower 
resolving powers of 10,000-30,000 are typically used 
to measure the large usini values of B-type stars (e.g. 
Halbedel 1996). We do not obtain usim estimates for 
these objects. 

For objects that do not have literature or synthetic 
spectrum fit RV estimates, we use the interstellar Na 
I D lines as an RV proxy. Finkenzeller & Jankovics 
(1984) were the first to show, using a sample of 10 
HAEBES, that the velocities of the molecular clouds in 
which HAEBES are embedded are very similar (mean 
residual of -3.1 ±3 km s“^) to the velocities of observed 
interstellar Ca H and Na I D velocities in the same ob¬ 
jects. They also showed that independently measured 
stellar RVs were similar to the molecular cloud veloci¬ 
ties, although with a larger mean residual of 8.5 ± 8 km 
s“^. Thus the Na I D velocities can, in general, provide 
a rough estimate of the stellar RV to within ^10 km s“^. 
As the mass flow diagnostics we are interested in are typ¬ 
ically several hundred km s“^ wide, this is adequate for 
our purposes. 

For objects in our sample for which we determine RVs 
based on synthetic spectrum fits, we compare these mea¬ 
surements to the RV determined from the object’s in¬ 
terstellar Na I D lines. RVs determined by the synthetic 
spectral fits from Alecian et al. (2013) are also used. The 
interstellar Na I RV estimate is made by manually point¬ 
ing at the approximate center of the interstellar absorp¬ 
tion. This is appropriate since the estimates are only 
accurate to ~8 km s“^ . An example of the Na I inter¬ 
stellar medium (ISM) absorption is shown in Figure 1 
for HD 142666. The Na I ISM absorption is very strong 
and narrow and is clearly distinguished from the photo¬ 
spheric profile. This is a particularly good example of 
the strength of using the Na I value as an estimate for 
the system RV: the Na I and synthetic model values dif¬ 
fer by only ~2.5 km s“^ . Further examples of the Na I 
ISM features can be seen in the individual object plots 
in the Appendix. 

The result of the comparison is plotted in Figure 2. 
Most objects show a velocity discrepancy of <10 km s“^ 

. The four objects that lie outside of the 20 km s“^ 
boundaries are all confirmed binaries and the large ve¬ 
locity differences may be due in part to their orbital mo¬ 
tion. We find a mean velocity difference of 8.1±3.6 km 
s“^ for the objects in Figure 2, where the uncertainty is 
the 95% confidence interval given by the t-distribution 
(Feigelson & Babu 2012). The mean difference drops to 
7.9 km s“^ if the four outliers are excluded, a negligible 
difference. The standard deviation of the absolute value 
of the difference between the interstellar RV values and 
those determined by the model fits is 8.5 km s“^ . We 
take this to be the uncertainty associated with the values 
determined from the Na I lines. 

Our sample thus confirms the findings of Finkenzeller 
& Jankovics (1984). For future HAEBE studies we sug- 
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gest that the interstellar Na I D line velocities be used 
for objects that, (1) do not allow for reliable photospheric 
fits either due to lack of absorption features or very high 
vsini, or (2) only ~10 km s“^ accuracy is needed in the 
radial velocity. All other physical information about the 


objects examined in this study (e.g. M*, i?*, M) are 
taken from the literature. The relevant references are 
given in Table 2. 


Table 2 

HAEBE optical sample physical parameters' 


Object ID 
(1) 

Spectral Type 
(2) 

^rad 

(km s“^) 

(3) 

M* 

(M©) 

(4) 

R* 

{Rq) 

(5) 

^’sin^ 
(km s“^) 
(6) 

log(M) 
(M© yr-i) 
(7) 

Disk 

Detected^ 

(8) 

n 

(9) 

References 

(10) 

AB Aur 

AO 

24.7 

2.50 

2.62 

116 

-6.85 

Y 

40 

1,2 

AE Lep 

B6 

28.0 







39 

BD+61 154 

B8 

-16.0 

3.40 

2.42 

112 


Y 

70 

1,6,9 

BF Ori 

A2 

22.0 

2.58 

3.26 

39 

<-8.00 

Y 


1,2,6,8 

BH Cep 

F5 

-5.6 

1.70 

2.40 

98 

<-8.30 

Y 

84 

2,38,40 

CQ Tau 

F2 

35.7 

2.93 

5.10 

98 

<-8.30 

Y 

29 

1,2,12,37 

DW CMa 

B3 






Y 


4,13 

CSC 04794-00827 

A? 

25.0 







43 

HD 141569 

AO 

35.7 

2.33 

1.94 

228 

-6.90 

Y 

55 

1,2,14 

HD 142666 

A5 

-9.7 

2.15 

2.82 

65 

-7.22 

Y 


1,3,15 

HD 163296 

A1 

-9.0 

2.23 

2.28 

129 

-7.16 

Y 


1,3,6,18 

HD 169142 

A7 

-0.4 

1.69 

1.64 

48 


Y 

0 

1,46 

HD 190073 

A1 

0.2 

2.85 

3.60 

4 


Y 

45 

1,19 

HD 203024 

A1 

-14.0 

2.80 

3.40 

162 


Y 

74 

1,44 

HD 244314 

A1 

22.5 

2.33 

2.07 

52 

-6.90 



1,3 

HD 244604 

A4 

26.8 

2.66 

3.69 

98 

-7.20 

Y 


1,3,21 

HD 245185 

A1 

16.0 

2.19 

1.85 

118 

-7.20 

Y 

81 

1,3,6,44 

HD 249879 

A2 

11.0 

4.00 

5.90 

249 

-8.00 



1,3 

HD 250550 

B8 

-22.0 

4.80 

3.50 

79 

-7.80 

Y 


1,3,13 

HD 287823 

AO 

-0.3 

2.50 

2.60 

10 


? 


1 

HD 290409 

A2 

80.3 

2.04 

1.75 

250 




1 

HD 290500 

A2 

29.0 

1.96 

1.68 

85 




1 

HD 290764 

A8 

26.3 

1.86 

2.30 

60 


Y 

32 

41 

HD 290770 

B9 

24.6 

2.86 

2.49 

240 




1 

HD 34282 

A3 

16.2 

1.59 

1.66 

105 

<-8.30 

Y 


1,2,18,22 

HD 35187 

A2 

27.0 

1.93 

1.58 

93 

-7.60 

Y 


1,3,24 

HD 35929 

FI 

21.1 

4.13 

8.10 

61.8 


Y 

56 

1,44 

HD 36408 

B8 

15.0 

4.10 

3.50 


<-8.00 

? 


2 

HD 37357 

A1 

21.4 

2.47 

2.83 

124 


Y 


1,47 

HD 37411 

B9 

14.1 





Y 

63 

33,44 

HD 38120 

B9 

28.0 

2.49 

1.91 

97 

-6.90 

Y? 

<8 

1,3,35 

HD 50083 

B4 

-0.5 

12.10 

10.0 

233 


7 


1 

HK Ori 

A3 

22.9 

3.00 

4.10 

20 

-5.24 

Y 


2,6 

IL Cep 

B4 

-39.0 



179 


? 


1 

IP Per 

A3 

13.7 

1.86 

2.10 

80 


7 


1 

IRAS 05044-0325 

B3 

25.2 







43 

IRAS 06071-1-2925 

B9 

12.5 







43 

IRAS 07061-0414 

BO 

43.0 







42 

IRAS 17481-1415 

AO 

-5.1 

2.50 






10 

IRAS 18306-0500 

B2 

33.1 

12.00 

6.80 



Y 

87 

40,44 

IRAS 18454-10250 

B1 

'^0 

19.00 

8.50 





42,13 

IRAS 19343-12926 

B1 

-10.7 







43 

LkHct 134 

B3 

-13.0 

6.00 

5.05 



Y 

58 

40,44 

LkHa 208 

A4 

13.6 

2.40 

3.20 



Y 

56 

40,44 

LkHo 233 

A3 

-9.10 

3.20 

4.45 



Y 

18 

40,44 

LkHo 257 

B5 

-15.2 

3.61 

2.68 



Y 

57 

44 

LkHo 324 

B8 

-10.5 

5.10 

6.30 



Y 

42 

40,44 

LkHo 339 

A1 

26.4 

3.18 

3.92 



Y 

40 

40,44 

MWC 1080 

B1 


17.4 

7.3 



Y 

83 

1,6,10 

MWC 120 

B9 

47.0 

3.94 

4.60 

120 

-6.85 

Y 


1,3 

MWC 137 

B1 






Y 

80 

6,10,13 

MWC 147 

B2 

18.6 

6.60 

4.80 

83 

-6.12 

Y 


1,3,13,48 

MWC 300 

B1 

-4.5 

7.37 

3.33 



Y 

62 

13,44 

MWC 361 

B4 

-7.0 

10.70 

10.40 

26 


Y 

26 

1,44 

MWC 480 

A4 

12.9 

1.93 

1.93 

98 

<-7.23 

Y 

37 

1,2,12 

MWC 610 

B3 

14.0 

8.00 

4.70 

219 


7 


1 

MWC 614 

AO 

15.1 



69 

-6.59 

Y 


1,3,28 

MWC 758 

A5 

17.8 

2.90 

4.40 

54 

-6.05 

Y 

21 

1,3,12 

MWC 778 

B2 

10.3 

10.00 

3.50 



Y 

75 

13,42,49 

MWC 863 

A1 

-5.0 

2.56 

2.89 

108 

-6.12 

Y 

38 

1,2,29 

NZ Ser 

B3 

-17.0 

10.10 

5.80 



Y 

5 

10,13 

R Mon 

B8 

22.3 

5.10 

12.00 



Y 

60 

10,13 

RR Tau 

AO 

15.3 

5.80 

9.30 

225 

-6.86 

Y 

44 

2,32,44 

SV Cep 

A1 

24.6 

2.62 

3.00 

180 


Y 

49 

1,44 

T Ori 

A3 

56.1 

3.13 

4.47 

147 

-6.60 

Y 


1,2,30 

UX Ori 

A1 

12.0 

6.72 

12.1 

221 

-7.18 

Y 

<8 

1,3,6,35 
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Table 2 — Continued 


Object ID 
(1) 

Spectral Type 
(2) 

'^rad 

(km s“^) 

(3) 

M, 

(Mq) 

(4) 

R* 

(^ 0 ) 

(5) 

vslni 
(km s“^) 
(6) 

log(M) 

(Mo yr-l) 
(7) 

Disk 

Detected*’ 

(8) 

n 

(9) 

References 

(10) 

UY Ori 

B9 

19.8 







4 

V1185 Tau 

A2 

20.0 



250 





V1578 Cyg 

Al 

-3.0 

5.90 

9.70 

199 


Y 


1,6 

V1685 Cyg 

B4 

-16.0 





Y 

41? 

30,34 

V1686 Cyg 

F8 

-12.2 

3.50 

14.00 



Y 


6,40 

V1787 Ori 

A5 

26.4 

1.88 

1.84 


-7.43 



45 

V1818 Ori 

B7 

26.4 



46 




43 

V346 Ori 

A7 

20.0 

1.72 

1.96 

116 

-6.90 

N? 


1,3,35 

V351 Ori 

A6 

15.0 

2.88 

4.38 

100 


7 


1 

V361 Cep 

B4 

-20.0 

8.11 

6.70 

278 


N? 


1,6 

V373 Cep 

B8 

-18.6 

5.30 

14.00 



Y 

45 

39,40 

V374 Cep 

B5 

-13.6 







51 

V380 Ori 

B9 

27.5 

2.87 

3.00 

7 

-5.60 

Y 


1,3,6 

V586 Ori 

Al 

31.0 

2.28 

1.94 

200 


Y 

81 

1,44 

V590 Mon 

B7 

17.7 

4.71 

5.86 




84 

40,44 

V599 Ori 

FO 

7.0 

1.83 

5.29 

96 


Y 

42 

43,44 

V718 Sco 

A4 

-3.6 

1.93 

2.25 

113 


Y? 

32? 

1,35 

V791 Mon 

B5 

-2.6 





7 


4 

VV Ser 

B7 

-2.0 

4.00 

3.10 

124 

-7.50 

Y 

46 

1,3,44 

WW Vul 

A2 

-4.0 

3.70 

5.40 

196 

-6.38 



1,2 

XY Per 

A2 

2.0 

1.95 

1.65 

224 

-7.02 

7 


1,3 

Z CMa 

B9 

-27.0 

3.80 

3.20 


-6.72 

Y 


3,31 


Note. — References: l=Alecian et al. (2013), 2=Mendigutfa 
et al. (2011), 3=Donehew &; Brittain (2011), 4=Vieira et al. 
(2003), 5=Carmona et al. (2010), 6=Hillenbrand et al. (1992), 
7=Tang et al. (2012), 8=Finkenzeller Mundt (1984), 9=Bossier 
et al. (2011), 10=Alonso-Albi et al. (2009), ll=Guilloteau et al. 
(2011), 12=Gmlloteau et al. (2013), 13=Verhoeff et al. (2012), 
14=Thi et al. (2014), 15=Schegerer et al. (2013), 16=Eisner et 
al. (2009), 17=Chen et al. (2012), 18=Marinas et al. (2011), 
19=Ragland et al. (2012), 20=Okamoto et al. (2009), 21=Vink 
et al. (2002), 22=Natta et al. (2004), 23=Acke & van den An- 
cker (2004), 24=Oudmaijer et al. (1992), 25=Fuente et al. (2002), 
26=Millan-Gabet et al. (2001), 27=Perrin et al. (2006), 28=Liu 
et al. (2007), 29=Fukagawa et al. (2003), 30=Eisner et al. (2004), 
31=Schutz et al. (2005), 32=Garcia Lopez et al. (2006), 33=Mal- 
fait et al. (1998), 34=Hernandez et al. (2004), 35=Dent et al. 
(2005), 36=Corporon &; Lagrange (1999), 37=Mendigutia et al. 
(2011), 38=Liu et al. (2011), 39=Lee Chen (2009), 40=Manoj 
et al. (2006), 41=Kraus et al. (2013), 42=Reed (2003), 43=Sar- 
tori et al. (2010), 44=Liu et al. (2011), 45=Caratti o Garatti et al. 
(2012), 46=Osorio et al. (2014), 47=Juhasz et al. (2010), 48=Bag- 
noli et al. (2010), 49=Perrin et al. (2009), 50=Acke et al. (2005), 
51=The, de Winter, & Perez (1994) 

^ All spectral types, values of M*, R*, M, i, and disk detections 
are taken from the literature. Values of the radial velocity and 
usim determined in this study are given in bold. Typical radial 
velocity uncertainties are ~5-10 km s“^; typical usim uncertain¬ 
ties are ~5-10 km s~^ although large values of usini can have 
uncertainties as large as 30-100 km s~^. 

^ Disk detections are indicated to provide context for the line 
profiles presented in this study. Question marks (?) indicate that 
the measurements are uncertain as to whether or not a disk is 
present. The detections are taken from the literature and are 
based mainly on millimeter and near-IR flux measurements (see 
given references). 



Velocity (km s"‘) 


Figure 1. Comparison of the interstellar Na I RV determination 
with the value calculated from fitting a synthetic spectrum to pho- 
tospheric absorption lines for HD 142666. The observed profile is 
shown in black; the broadened spectroscopic standard HR 3221 is 
over plotted in red. The RV value estimated from the ISM feature 
is shown with the vertical red line. The photospheric model fit RV 
is shown with the vertical blue line. The ±8.5 km s~^ uncertainty 
for the Na I value is shown with the vertical green dashed lines. It 
is clear that manually selecting the center of the ISM absorption re¬ 
sults in negligible errors compared with the systematic uncertainty 
of using the Na I to estimate the system RV. 


4. PROFILE CLASSIFICATION 

Line profiles are classified into six morphological 
groups: P-Cygni (PC), inverse P-Cygni (IPC), double- 
peaked emission (DP), single-peaked emission (E), ab¬ 
sorption (A), and featureless (F). For the purpose of 
generating mass flow statistics, pure absorption profiles 

















Figure 2. Radial velocities calculated from synthetic model pho¬ 
tosphere fits versus Na I D interstellar absorption velocities. The 
solid line represents equal velocities. The dashed lines moving out¬ 
ward from the solid line represent differences between the velocities 
of 10 and 20 km s“^ . The mean deviation is 8.1 km s~^ . We note 
that the four objects (HD 50083, IL Cep, MWC 614, and SV Cep) 
lying outside of the 20 km s“^ boundary are all known binaries. 
The large velocity discrepancies in these objects may be due to 
their orbital motion. 

(A) are further grouped as red or blue absorption based 
on the maximum absorption velocity: Vred > 0 km s“^ 
and Vbiue < 0 km s“^. The objects counted as showing 
blue(red)-shifted absorption either have a PC(IPC) pro¬ 
file or an A profile that is blue(red)-shifted. Objects are 
classified as E or DP if the profile is in emission (single 
and double peaked, respectively) and there is no clear ev¬ 
idence of absorption below the local continuum. Profiles 
are labeled as F if they show no detectable departure 
from the estimated photospheric line proHle. We note 
that a small number of the PC profiles, most of them 
Ha profiles, show blue-shifted absorption that does not 
extend below the stellar continuum, i.e., the absorption 
is superposed on the wing of the emission prohle (see, for 
example, the Ha profile in Figure 3). While these are not 
strictly P-Cygni profiles, we have chosen to extend the 
PC classification to these objects since it is the closest 
morphology group that describes the profile. 

In order to accurately classify the morphologies, each 
profile was compared with a rotationally broadened spec¬ 
troscopic standard of a similar effective temperature. 
Objects with unknown risim values did not have their 
comparison spectra broadened. For most objects, this 
allowed weak emission and absorption features to be 
identified in lines with significant photospheric compo¬ 
nents. This comparison is most important for the Balmer 
lines since most of our objects have strong photospheric 
Balmer line absorption. An example of this effect is 
shown in the top panel of Figure 5 where the core of 
the Hy profile of HD 190073 is clearly filled in by emis¬ 
sion. A second example is shown in the bottom panel 
of Figure 5 for a more complex line morphology at H/3 
for HD 290770. In this case the blue-shifted absorption 
at ~-100—200 km s“^ can clearly be seen superimposed 


on the broad emission profile. The comparison with the 
standard makes it obvious that this is non-photospheric 
absorption. To be clear, we note that the profiles shown 
in Figure 3 and Figure 4 show only the continuum nor¬ 
malized profiles-no photospheric subtraction was done 
to make these plots even though photospheric compar¬ 
isons were made while classifying the profiles. The pho¬ 
tospheric comparisons are done by eye and, in general, 
it is obvious when a non-photospheric contribution to 
the line profile is present. Thus no quantitative compari¬ 
son of the photospheric and circumstellar components is 
performed. The approximate comparisons prevent gross 
misclassification of line profiles, i.e. identifying circum¬ 
stellar absorption when the absorption actually forms in 
the photosphere. Thus objects will often show line pro¬ 
files that appear to show structure but are in fact due 
entirely to photospheric absorption. Examples of an AO 
and B3 spectroscopic standard are shown in Figure 6 and 
Figure 7. 

In some objects there is significant contamination of 
the He I 3889 A, Mg I 5167 A, Ca H 3969 A, and Ca 
H 8662 A spectral lines by neighboring features. This 
is especially severe at Mg I 5167 A (e.g. Z CMa and 
DW CMa), where we tentatively identify the nearby Fe I 
5168.89 A line (-1-81.8 km s“^), and Ca H 3969 A and He 
I 3889 A which are often overwhelmed by the nearby He 
and HC lines, respectively. For example, in Figure 3 the 
He I 3889 A line appears to show a PC morphology. How¬ 
ever, this profile more likely due to the nearby H^ line. 
Thus this line is listed as “contaminated” since the true 
profile is obscured by the H^ line. The N H 6548 A line 
lies in the far blue wing of the Ha and often appears to 
show structure when in fact the line emission and shape 
are due entirely to Ha (e.g.. Figure 3). If there are other 
lines of the same element that show a non-featureless 
line profile, e.g. the Mg I 5173 or 5184 A lines in the 
case of Mg I 5167 A, the line is categorized as ’C’ for 
’contaminated’. If the other lines of the same element 
are featureless, we assume a featureless spectrum for the 
feature in question. The rest velocities of these contami¬ 
nating lines are marked with a dotted line in the individ¬ 
ual profiles plots. With the exception of Ca H 8662 A, 
the wavelength regions containing these diagnostics are 
mostly featureless in the whole sample. Contaminated 
profiles comprise 16% of the He I 3889 A profiles, 3% of 
the Mg I 5167 A profiles, 7% of the Ca H K profiles, and 
9% of the Ca H 8662 A profiles. Thus the contamination 
is not problematic concerning the classihcation of most 
profile morphologies. 

Examples of the extracted profiles and their catego¬ 
rization are shown in Figure 3 and Figure 4. Plots of 
the individual line profiles for all objects are given in 
Appendix A. The profile classihcation is shown in the 
upper right corner of each plot window; the spectral fea¬ 
ture being shown is listed in the upper left. Spectral lines 
for which there is no data are labeled as such in the plots. 
These lines either fell in between orders on the detector 
or the S/N was too low for analysis. Features with rest 
wavelengths below 5000 A are plotted using a 5 km s“^ 
(2 pixel) bin window in order to enhance the S/N. The 
bottom-right panel of each hgure shows a bar plot of the 
number of line prohles in each morphology category. We 
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note that only the PC and IPC bars are shown in color. 
The bar showing the count for pure red- or blue-shifted 
absorption (i.e., the A profiles) is shown in black even 
though the subdivision of the A profiles into red- and 
blue-shifted absorption is included in the statistics. 

5. MORPHOLOGY STATISTICS 

5.1. Pre-main sequence evolutionary groups 

HAEBES display a wide variety of SEDs indicating 
that, as a group, they represent objects at different 
phases of pre-main sequence evolution (Meeus et al. 
2001; Malfait et al. 1998; Hillenbrand et al. 1992). An 
accurate calculation of mass flow statistics would only 
include objects that are currently in the same pre-main 
sequence phase as one another. However, this is challeng¬ 
ing due to the relatively small number of HAEBES each 
specific evolutionary phase would contain. Studies of this 
sort will necessarily suffer from smaller sample sizes in 
each group thus increasing the uncertainties associated 
with the mass flow statistics. 

In order to eliminate any objects from our sample that 
no longer show evidence of inner disk material, and thus 
should not show evidence of mass accretion and mass 
loss, we employ the color-color classification of Hillen¬ 
brand et al. (1992). The J — H vs H — K classification 
is a proxy for the full SED analysis presented in Hil¬ 
lenbrand et al. (1992). Objects with J — H<0.3 and 
H — K<0.3, the Group HI objects from Hillenbrand et 
al. (1992), likely have insufficient material in their in¬ 
ner disks to be actively accreting and ejecting material. 
These objects should not be included in any comparisons 
to CTTSs since all CTTSs show strong evidence for ac¬ 
cretion and mass loss. 

We have collected J, H, and K data from the liter¬ 
ature for our observed sample in order to identify any 
potential Group HI objects. Optical extinctions are also 
collected for some objects. The J — H and H — K col¬ 
ors for our sample are plotted in Figure 8. Objects with 
de-reddened colors are plotted as green circles while ob¬ 
jects without known Ay values are plotted as red squares. 
Reddening vectors for the standard i?=3.1 extinction law 
are plotted as solid black lines. An approximation of the 
main sequence for dwarf spectral types B8 through M6 
is represented by the dashed-dotted line. Main sequence 
objects are required to remain within the boundaries of 
the reddening vectors. It can be seen that most of our 
objects lie well to the right of the main sequence region, 
confirming the existence of significant excess emission. 
The Group HI objects from Hillenbrand et al. (1992) all 
exist in region J — H<0.3, H — K<0.3. We use the 
box defined by the intersection of these boundaries as 
our criteria for identifying Group HI objects in our sam¬ 
ple. The box is plotted in Figure 8 using dashed lines. 
There are 9 objects (AE Lep, HD 141569, HD 203024, 
HD 35929, HD 36408, HD 50083, IL Cep, MWC 610, and 
V361 Cep) that lie within the Group HI box. These ob¬ 
jects are marked with a f in Table 4 and will be excluded 
from the analysis presented in the rest of the paper. This 
reduces the sample to 78 objects. 

5.2. Calculating the mass flow incidence 

In order to exclude objects that are not unambigu¬ 
ously experiencing some sort of outflow or accretion, we 


do not include pure emission lines (double or single- 
peaked) as evidence of mass flows. While emission lines 
in many species can form in a wind or accretion flow, 
they can also form in static geometries such as extended 
stellar and disk atmospheres or remnant circumstellar en¬ 
velopes. On the other hand, doppler shifted absorption 
is an unambiguous sign of mass motion. The CTTS com¬ 
parison statistics that are discussed in subsection 5.5 are 
calculated using the same criteria. Thus we only count 
red and blue-shifted absorption below the local contin¬ 
uum (e.g. the blue-shifted absorption seen at H/3 in the 
bottom panel of Figure 5) as indications of mass flows 
when calculating the final statistics in Table 3 and Ta¬ 
ble 4. For objects with multiple observations (^20% of 
the sample, see Table 1) we search for absorption signa¬ 
tures in all of the exposures. Thus the computed statis¬ 
tics include objects that, for example, show blue-shifted 
absorption in a line for one exposure but do not show the 
same morphology in a separate exposure. This enables 
us to account for variability, in a limited way, in at least 
some objects. 

The morphology statistics for the mass flow diagnostics 
are presented in Table 3 and a summary of each object’s 
profile classifications is given in Table 4. In Table 3 the 
line diagnostic with the highest occurrence of each profile 
type is highlighted in bold font. For each line the total 
number of object profiles examined is give in column 2. 
The number of each type of morphology is listed and the 
fraction of the total number of profiles is given in paren¬ 
theses. In column 8 of Table 4 the total number of line 
diagnostics examined is given for each object. The num¬ 
ber of contaminated profiles for each object is not shown 
although they are included in the total number of lines 
in column 8. It is immediately obvious from Table 3 that 
most of the diagnostics show featureless profiles. In fact, 
Ha, H/?, and Ga H A8498 are the only lines for which 
a majority of the sample displays some form of activity. 
The lines with the highest incidence of circumstellar fea¬ 
tures displaying PG and IPG morphologies are H/3 (30%) 
and He I A5876 A (14%), respectively. The Na I D lines 
show the highest rate (12%) of pure absorption profiles. 
After the Balmer lines, the next best spectral line for ob¬ 
serving a PC morphology is Fe H A4924 which displays 
PC characteristics in 11% of the object profiles. 

The last row in Table 4 displays the final mass flow 
statistics for our sample. The uncertainties are the 68% 
confidence intervals calculated using Wilson’s score test 
for large samples (Wilson 1927). These intervals are very 
close to the values obtained assuming Poisson statistics 
in the distribution of observed incidence. However, the 
Wilson score test results in a weighted average towards 
the center of the distribution. We note that the given 
incidence is not the adjusted fraction given by the score 
test and thus the confidence intervals may be asymmetric 
around the value. The occurrence of both inflow and 
outflow signatures in our sample are fairly low: only 26% 
of our objects have at least one profile displaying red- 
shifted absorption; only 37% have at least one profile 
displaying blue-shifted absorption. 

5.3. Comparison to previous HAEBE studies 

A comparison of our mass flow statistics to previous 
HAEBE studies is presented in Table 5. We find simi- 
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Velocity (km s ') 


Figure 3. Extracted line profiles for Z CMa. The continuum, which is normalized to 1.0, is marked with a horizontal dashed line. The 
stellar rest velocity is marked with a vertical dashed line. The profile classification is noted in the upper-right of each plot window. The 
plotted spectral line is listed in the upper—left. The bottom right-most plot window is a bar plot of the profile types for the object. The 
featureless profile count is excluded from the bar plot. Z CMa displays strong outflow signatures in numerous lines. The Ha vertical plot 
range is abbreviated to show the blue—shifted absorption. The [Nil] 6548 A appears to show emission but is actually showing the far blue 
wing of the Ha line. Although the Na I doublet lines should show scaled versions of the same profile morphology, the blue—shifted absorption 
from the 5896 A profile decreases the emission peak of the 5889 A profile to below the local continuum, resulting in a classification of ‘A’ 
for the 5889 A line. Note the vertical dotted lines denoting the rest wavelengths of lines that potentially contaminate the line of interest 
(e.g., H(^ in the He I 3889 A panel). 
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Figure 4. Extracted line profiles for MWC 120. The format is the same as Figure 3. Note the strong IPC morphology in the He I 5876 
A line. 
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Figure 5. A comparison of the H 7 profile of HD 190073 with 
a spectroscopic standard of similar Teff (top panel). The object 
profile is plotted in black and the standard photospheric profile is 
plotted in red. The core of the line is clearly filled in by emission, 
placing this profile in the E morphology group. A similar compar¬ 
ison is shown for HD 290770 (bottom panel) at H/3. The standard, 
HR4386, has been rotationally broadened to match the usini of HD 
290770. The blue-shifted absorption can be seen superimposed on 
the emission profile, resulting in a profile classification of PC. 


lar rates of blue-absorption in our sample compared to 
both Finkenzeller & Mundt (1984, FM84) and Hamann 
& Persson (1992, HP92). Both Vieira et al. (2003, V03) 
and Reipurth et al. (1996, R96) find rates of blue-shifted 
absorption that are significantly lower than ours (as well 
as FM84’s and F[P92’s). We attribute this large discrep¬ 
ancy to the fact that they only observe Hct which tends to 
mask sub-continuum absorption due to strong emission 
components from multiple sources (e.g. emission from 
accretion, stellar winds, nebular background). We note, 
however, that using solely Ha in our sample results in a 
blue-shifted absorption fraction of 26%l4, ^ value in be¬ 
tween the low rates of Vieira et al. (2003) and Reipurth 
et al. (1996) and our final fraction (37%). This may be 
due to our classification of profiles as P-Cygni if the blue 
feature is clearly absorption and not a second peak, even 
if the absorption does not extend below the local emis¬ 
sion continuum (e.g. Figure 5). Our results, and those 
of FM84 and HP92, show that in order to recover accu¬ 
rate outflow statistics diagnostics other than Ha should 
be examined. The line with the largest occurrence of 
blue-absorption in our sample is H,0 (30%), although 
the incidence in Ha is similar at 26%. In fact, all objects 
in our sample that show P-Cygni type proHles in lines 
other than Ha and H/3 also show P-Cygni morphologies 
in either Ha or H/3. Thus it appears H/3 is a very strong 
stand-alone diagnostic of outflows in HAEBES. 


Table 3 

Line morphology statistics 


Line ID 
( 1 ) 

Total 

( 2 ) 



Profile morphology 



PC 

(3) 

IPC 

(4) 

DP 

(5) 

E 

( 6 ) 

A 

(7) 

F 

( 8 ) 

Ha 

78 

20(.26) 

5(.06) 

34(.44) 

17(.22) 

o(.oo) 

2(.03) 

H/3 

74 

22(.30) 

4(.05) 

15(.20) 

18(.24) 

5(.07) 

10(.14) 

Hy 

61 

12 (. 20 ) 

0 {. 00 ) 

3(.05) 

14(.23) 

3(.05) 

29(.48) 

H<5 

51 

5(.10) 

O(.OO) 

1 (. 02 ) 

8(.16) 

5(.10) 

32(.63) 

He I A3889 

49 

O(.OO) 

0600) 

o(.oo) 

1 (. 02 ) 

O(.OO) 

40(.82) 

He I A5876 

78 

1 {. 01 ) 

11(.14) 

l(.Ol) 

11(.14) 

If.Ol) 

53(.68) 

He I A6678 

71 

l(.Ol) 

l(.Ol) 

o(.oo) 

2(.03) 

l(.Ol) 

66(.93) 

Mg I A5167 

71 

O(.OO) 

O(.OO) 

O(.OO) 

O(.OO) 

O(.OO) 

69(.97) 

Mg I A5173 

71 

0 {. 00 ) 

0600) 

0600) 

3604) 

Of.OO) 

68(.96) 

Mg I A5184 

71 

l(.Ol) 

O(.OO) 

O(.OO) 

5(.07) 

O(.OO) 

65(.92) 

Ca H A3933 

57 

2(.03) 

l(.Ol) 

0600) 

l(.Ol) 

4(.05) 

65(.89) 

Ca H A3969 

42 

O(.OO) 

O(.OO) 

O(.OO) 

3(.07) 

1 {. 02 ) 

35(.83) 

Ca H A4227 

71 

l(.Ol) 

O(.OO) 

O(.OO) 

o(.oo) 

O(.OO) 

70(.99) 

Ca H A8498 

27 

2(.07) 

1(.04) 

2(.07) 

9(.33) 

O(.OO) 

13(.48) 

Ca H A8662 

70 

2(.03) 

2(.03) 

3(.04) 

22(.31) 

O(.OO) 

35(.50) 

Na I A5889 

78 

4(.05) 

l(.Ol) 

O(.OO) 

19(.24) 

9(.12) 

45(.58) 

Na I A5896 

78 

4(.05) 

1 {. 01 ) 

O(.OO) 

18(.23) 

8 (. 10 ) 

47(.60) 
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Figure 6. Extracted line profiles for the spectroscopic standard HR 4781. 
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Figure 8. J — H and H — K colors for our sample. Redden¬ 
ing vectors are plotted as solid lines. The approximate location 
of the main sequence is indicated by the dashed—dotted line. The 
area within the dashed lines contains objects that are most likely 
evolved HAEBES without any significant IR excess. Objects with 
de-reddened colors are plotted as green circles; objects without 
reddening corrections are plotted as red squares. There are 9 ob¬ 
jects contained in the box that we exclude from further analysis. 

Table 3 — Continued 


Profile morphology 



Line ID 
(1) 

Total 

(2) 

PC 

(3) 

IPC 

(4) 

DP 

(5) 

E 

(6) 

A 

(7) 

F 

(8) 
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Profile morphology 




Line ID 

Total 

PC 

IPC 

DP 

E 

A 

F 



(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 
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Profile morphology 




Line ID 

Total 

PC 

IPC 

DP 

E 

A 

F 



(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 


N II 

A5755 

78 

0 {.00) 

O(.OO) 

O(.OO) 

3(.04) 

0 {.00) 

75(.96) 


N II 

A6548 

78 

0 ^00) 

O(.OO) 

0 ^00) 

2^03) 

0 ^00) 

76^97) 


N II 

A6583 

40 

0 ^00) 

O(.OO) 

0 ^00) 

6(.15) 

0 ^00) 

34^85) 


O I 

A5577 

78 

0 ^00) 

0 ^00) 

O(.OO) 

6(.08) 

0 ^00) 

72^92) 


O I 

A6300 

77 

0 ^00) 

0 ^00) 

0 ^00) 

43(.56) 

0 ^00) 

34^44) 


O I 

A6363 

64 

0 ^00) 

0 ^00) 

0 ^00) 

20(.31) 

0 ^00) 

44^69) 


O III] A5007 

71 

0 ^00) 

0 ^00) 

O(.OO) 

3(.04) 

0 ^00) 

68^96) 

Fe I A5498 

78 

0 ^00) 

0 ^00) 

0 ^00) 

l(.Ol) 

0 ^00) 

77^99) 

Fe II A4924 

64 

7{.ll) 

2^03) 

1 (.02) 

7(.ll) 

3^05) 

44^69) 

Fe II A4352 

55 

3(.05) 

0 ^00) 

O(.OO) 

3(.05) 

0 ^00) 

49^89) 


Table 4 

Object profile classification summary 


Object ID 
(1) 

Npc 

(2) 

Nipc 

(3) 

Nop 

(4) 

Ne 

(5) 

Na 

(6) 

Nf 

(7) 

nL 

(8) 

Red abs. 
(9) 

Blue abs. 
(10) 

AB Aur 

4 

0 

0 

5 

1 

16 

26 

No 

Yes 
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Table 4 — Continued 


Object ID 
(1) 

(2) 

Nipc 

(3) 

Nop 

(4) 

Ne 

(5) 

Na 

(6) 

Np 

(7) 

(8) 

Red abs. 
(9) 

Blue abs. 
(10) 

AE Lept 

0 

0 

0 

5 

0 

7 

12 

No 

No 

BD+61 154 

6 

0 

0 

9 

0 

10 

26 

No 

Yes 

BF Ori 

0 

2 

1 

0 

8 

15 

27 

Yes 

No 

BH Cep 

1 

0 

0 

0 

3 

22 

26 

No 

Yes 

CQ Tau 

0 

0 

1 

2 

0 

24 

27 

No 

No 

DW CMa 

5 

0 

0 

7 

0 

12 

25 

No 

Yes 

CSC 04794-08227 

0 

0 

0 

6 

0 

20 

26 

No 

No 

HD 141569t 

0 

0 

2 

1 

0 

27 

30 

No 

No 

HD 142666 

0 

1 

1 

1 

0 

25 

28 

Yes 

No 

HD 163296 

1 

1 

0 

3 

1 

21 

27 

Yes 

Yes 

HD 169142 

0 

0 

0 

2 

0 

23 

25 

No 

No 

HD 190073 

2 

0 

0 

12 

0 

14 

29 

No 

Yes 

HD 203024t 

0 

0 

0 

1 

0 

26 

27 

No 

No 

HD 244314 

2 

0 

0 

7 

0 

20 

29 

No 

Yes 

HD 244604 

2 

0 

0 

5 

0 

19 

26 

No 

Yes 

HD 245185 

0 

0 

1 

3 

2 

20 

26 

Yes 

No 

HD 249879 

2 

0 

0 

4 

0 

18 

24 

No 

Yes 

HD 250550 

2 

0 

0 

4 

0 

6 

12 

No 

Yes 

HD 287823 

1 

0 

0 

1 

0 

24 

26 

No 

Yes 

HD 290409 

0 

0 

1 

1 

0 

24 

26 

No 

No 

HD 290500 

0 

0 

1 

0 

0 

25 

26 

No 

No 

HD 290764 

0 

0 

0 

2 

0 

24 

26 

No 

No 

HD 290770 

1 

0 

0 

5 

0 

20 

26 

No 

Yes 

HD 34282 

0 

0 

1 

0 

0 

25 

26 

No 

No 

HD 35187 

0 

2 

0 

1 

0 

24 

27 

Yes 

No 

HD 35929t 

0 

0 

1 

1 

0 

24 

26 

No 

No 

HD 36408t 

1 

0 

0 

5 

0 

19 

25 

No 

Yes 

HD 37357 

2 

0 

0 

1 

0 

23 

26 

No 

Yes 

HD 37411 

0 

0 

1 

1 

0 

10 

12 

No 

No 

HD 38120 

0 

1 

1 

2 

0 

8 

12 

Yes 

No 

HD 50083t 

0 

0 

3 

1 

0 

21 

26 

No 

No 

HK Ori 

0 

0 

2 

3 

0 

21 

26 

No 

No 

IL Cept 

0 

0 

2 

0 

0 

24 

26 

No 

No 

IP Per 

0 

0 

0 

4 

0 

22 

26 

No 

No 

IRAS 05044-0325 

0 

0 

2 

1 

0 

22 

25 

No 

No 

IRAS 06071-1-2925 

2 

0 

0 

5 

0 

18 

26 

No 

Yes 

IRAS 07061-0414 

0 

0 

0 

0 

0 

25 

25 

No 

No 

IRAS 17481-1415 

0 

0 

0 

2 

0 

18 

21 

No 

No 

IRAS 18306-0500 

0 

0 

0 

0 

0 

23 

23 

No 

No 

IRAS 18454-1-0250 

0 

0 

0 

1 

0 

19 

20 

No 

No 

LkHa 134 

2 

0 

0 

6 

0 

18 

26 

No 

Yes 

LkHa 208 

0 

0 

1 

2 

0 

24 

27 

No 

No 

LkHa 233 

2 

0 

0 

4 

0 

17 

23 

No 

Yes 

LkHa 257 

0 

0 

2 

0 

0 

20 

22 

No 

No 

LkHa 324 

0 

0 

2 

0 

0 

21 

23 

No 

No 

LkHa 339 

0 

0 

0 

3 

0 

19 

22 

No 

No 

MWC 1080 

9 

0 

0 

4 

0 

10 

25 

No 

Yes 

MWC 120 

0 

3 

3 

7 

0 

12 

26 

Yes 

No 

MWC 137 

0 

0 

0 

13 

0 

10 

24 

No 

No 

MWC 147 

0 

0 

3 

3 

0 

16 

22 

No 

No 

MWC 300 

0 

0 

2 

12 

0 

12 

26 

No 

No 

MWC 361 

0 

0 

3 

1 

0 

21 

26 

No 

No 

MWC 480 

6 

0 

1 

1 

3 

13 

25 

No 

Yes 

MWC 610+ 

0 

0 

2 

0 

0 

27 

29 

No 

No 

MWC 614 

0 

0 

0 

2 

0 

24 

26 

No 

No 

MWC 758 

2 

0 

0 

4 

0 

20 

26 

No 

Yes 

MWC 778 

0 

0 

2 

3 

0 

8 

13 

No 

No 

MWC 863 

2 

0 

0 

3 

0 

20 

25 

No 

Yes 

NZ Ser 

0 

0 

0 

7 

0 

13 

21 

No 

No 

R Mon 

1 

0 

1 

7 

0 

13 

22 

No 

Yes 

RR Tau 

0 

1 

2 

3 

0 

21 

27 

Yes 

No 

SV Cep 

0 

0 

1 

1 

0 

23 

25 

No 

No 

T Ori 

0 

1 

2 

5 

0 

19 

27 

Yes 

No 

UX Ori 

0 

5 

0 

3 

0 

15 

23 

Yes 

No 

UY Ori 

0 

0 

1 

2 

2 

7 

12 

No 

Yes 

VI185 Tau 

1 

0 

0 

2 

0 

24 

27 

No 

Yes 

V1578 Cyg 

7 

0 

0 

2 

2 

14 

26 

No 

Yes 

V1685 Cyg 

1 

1 

4 

3 

0 

20 

29 

Yes 

Yes 

V1686 Cyg 

0 

0 

2 

2 

0 

16 

19 

No 

No 

V1787 Ori 

0 

0 

1 

1 

0 

21 

23 

No 

No 

V1818 Ori 

0 

0 

5 

1 

2 

16 

24 

Yes 

No 

V346 Ori 

0 

1 

0 

0 

1 

24 

26 

Yes 

No 

V351 Ori 

0 

2 

0 

0 

2 

22 

26 

Yes 

No 

V361 Cept 

0 

0 

1 

3 

0 

21 

26 

No 

No 
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Table 4 — Continued 


Object ID 
(1) 

Tco 

(2) 

NiPC 

(3) 

Ndp 

(4) 

Ne 

(5) 

Na 

(6) 

Np 

(7) 

NLt 

(8) 

Red abs. 
(9) 

Blue abs. 
(10) 

V373 Cep 

3 

0 

0 

3 

0 

15 

22 

No 

Yes 

V374 Cep 

3 

0 

1 

0 

2 

22 

29 

No 

Yes 

V380 Ori 

0 

0 

0 

8 

0 

5 

14 

No 

No 

V586 Ori 

0 

1 

1 

6 

2 

16 

26 

Yes 

No 

V590 Mon 

0 

0 

2 

4 

0 

20 

26 

No 

No 

V599 Ori 

0 

1 

0 

2 

0 

15 

18 

Yes 

No 

V718 Sco 

0 

1 

0 

0 

1 

21 

23 

Yes 

No 

V791 Mon 

2 

0 

0 

6 

0 

5 

13 

No 

Yes 

VV Ser 

0 

2 

3 

3 

0 

17 

25 

Yes 

No 

WW Vul 

0 

2 

1 

3 

2 

20 

28 

Yes 

No 

XY Per 

0 

1 

1 

0 

5 

18 

26 

Yes 

No 

Z CMa 

13 

0 

0 

5 

1 

7 

27 

No 

Yes 

Total fraction"^: 

0.36 

0.23 

0.45 

0.83 

0.22 



U.^d_q Q5 

0 Q7+0.06 


^ Total number of line profiles that show the specified morphology. 

^ Total number of line profiles examined for each object. 

^ Total fraction of objects that display the specified profile or absorp¬ 
tion type. 


Table 5 

Comparison to previous studies 


Study 

(1) 

Ntot 

(2) 

(3) 

% RA 
(4) 

Nr A 

(5) 

% BA 
(6) 

Present work 

78 

20 

26ti 

29 

37l® 

FM84 

46 

2 

4t^ 

16 


HP92 

32 

2 

6 l« 

11 


R96 

18 

1 

6 t« 

2 

lltl 

VOS 

131 



18 



The difference in red-absorption statistics between our 
study and the others listed in Table 5 is more pro¬ 
nounced. We find that 26% of our objects show red- 
shifted absorption in at least one line while the studies of 
FM84, HP92, and R96 show roughly a 4-6% occurrence. 
The discrepancy is easily explained for our study versus 
those of FM84 and R96: the occurrence of red-shifted 
absorption at Ha in our study is ~6%, comparable to 
that found by FM84 and R96 both of which primarily 
examine Ha. FM84 also examine the Na I D lines but 
both their sample and ours show a negligible occurrence 
of red-shifted absorption in those lines. 

The discrepancy with HP92 is harder to explain. A 
large percentage (15%) of our objects show red-shifted 
absorption at He I A5876. The rest of the red-shifted 
absorption detections result mainly from Ha, H/3, and 
Ca H IRT lines. HP92 also examine these lines but they 
detect red-shifted absorption in only 2 out of 32 objects, 
a detection rate of ~6%. This might be explained by 
the fact that HP92 only observe 6 objects at He I A5876 
and thus, to obtain a detection rate at A5876 similar to 
that of our sample, they would, on average, detect red- 
shifted absorption in ~1 object, which they do. The 
non-detections in the other lines is more puzzling but 
may be accounted for by the lower resolution of their ob¬ 
servations (11-85 km s“^). These medium-resolution ob¬ 
servations may mask weak red-shifted absorption signa¬ 
tures. We note that the red-shifted absorptions observed 


at Ha and H/3 in our sample tend to appear as deeper 
absorptions superimposed on a photospheric absorption 
core. However, when these lines are carefully compared 
to a rotationally broadened standard it is clear that the 
red-shifted absorptions are non-photospheric and that 
the profile morphology is inverse P-Cygni. It is not clear 
that HP92 compare their objects to photospheric stan¬ 
dards. If not, this may also contribute to the discrepancy. 

By comparison to the studies listed above, our anal¬ 
ysis demonstrates the necessity of large sample, high- 
spectral resolution studies of outflow and inflow diagnos¬ 
tics in HAEBES in order to gather accurate statistics 
of mass flows around these objects. The H/3 diagnos¬ 
tic shows the largest percentage of P-Cygni profiles, al¬ 
though other lines (e.g. Ha, Hy, Fe I 4924) also show 
non-negligible P-Cygni detection fractions. The He I 
A5876 line is by far the most sensitive line for detecting 
inverse P-Cygni morphologies. This is not surprising 
since the He I A10830 line is also a very efficient diagnos¬ 
tic for detecting mass Hows (Edwards et al. 2006; Cauley 
& Johns-Krull 2014) and the A5876 transition lies im¬ 
mediately above the upper energy level of the A10830 
transition, preceding A10830 in a recombination cascade 
sequence. The Na I D lines are also relatively sensitive to 
outflows in HAEBES, although the lines more frequently 
display pure emission profiles. 

5.4. Comparison to He I X10830 data 

The line profile statistics (Table 6) for both inflows and 
outflows from our optical data set are similar to those 
from the He I A10830 data set, although the fractions 
are higher in the He I A10830 sample. 

We use contingency tests to compare the fractions of 
red and blue-shifted absorption between the two studies. 
A contingency test of the measured fraction of a charac¬ 
teristic in two groups gives the probability that the null 
hypothesis is true, i.e. that the intrinsic fraction is the 
same in both groups (Feigelson & Babu 2012). This prob¬ 
ability is referred to here as the p-value. A high p-value 
indicates a high probability of having obtained the dif¬ 
ference in the measured fractions by chance, i.e. the null 
hypothesis cannot be ruled out; a low p-value indicates a 
low probability that the difference in the measured frac¬ 
tions was obtained accidentally, i.e. the difference is real 
and the null hypothesis can probably be ruled out. A 
contingency test comparing the He I A10830 sample to 
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the optical sample shows that neither the blue-shifted 
(p =0.85) nor the red-shifted (p =~1.0) absorption inci¬ 
dences are signihcantly different. 


Table 6 

Comparison to He I A10830 HAEBE statistics 



Optical data 

He I A10830 data 


(1) 

(2) 

Ntot 

78 

48 

Red—shifted abs. 

26%1® 

27%t6 

Blue—shifted abs. 

37%tl 

40%j;/ 


Although the overall morphology statistics are simi¬ 
lar between the optical and He I A10830 data, the He I 
A10830 line is much more efficient at tracing mass flows 
around HAEBES than any single optical diagnostic. In 
addition, He I A10830 is generally less contaminated by 
sources other than nearby circumstellar material due to 
the high excitation temperatures needed to populate the 
lower energy level (~20,000 K; Dupree et al. 2005). This 
is not the case for other popular diagnostics such as Ha, 
H/3, and the Na I D lines which can frequently display 
strong emission signatures from chromospheres, nebular 
background, and interstellar absorption in the case of Na 
I D. The tendency of the He I A10830 line to remain in 
its lower energy level until a photon is absorbed accounts 
for the large opacity needed to generate the blue-shifted 
absorption commonly seen in the data from Paper I. 

Figure 9 compares the maximum absorption velocities 
of He I A10830 and Ha or H/3 for objects that display 
blue-shifted absorption in both He I A10830 and Ha (red 
diamonds) and H/3 (blue crosses). This figure shows that 
He I A10830 traces approximately the same outflow ve¬ 
locities as Ha and H,5, although the velocities in H/3 are 
generally smaller when compared to He I A10830. This 
could be due to the larger optical depth of He I A10830 
and Ha relative to H/3, allowing photons to be absorbed 
up to larger fractions of the terminal wind velocity. The 
number of objects with both optical and He I A10830 
data that also show blue-shifted absorption is small and 
thus it is not possible to make a strong case for choos¬ 
ing one of the optical diagnostics over He I A10830 based 
on a velocity comparison. Each line appears to trace 
similar material in the wind. Thus for future studies 
concerned with detecting the presence of outflows and 
accretion around HAEBES, He I A10830 could serve as 
a primary diagnostic due to its high efficiency in detect¬ 
ing mass flows. In addition, it is relatively free of pho- 
tospheric and interstellar contaminants making the cir¬ 
cumstellar nature of the line morphologies unambiguous. 
In terms of mass flow statistics and outflow velocities, 
it appears that optical data does not provide any addi¬ 
tional information over that provided by the He I A10830 
line. 

5.5. Comparison to CTTSs 

To our knowledge, the study performed by Alencar & 
Basri (2000, hereafter ABOO) is the only current high- 
resolution, moderate sample size investigation of CTTSs 
that examines many different spectral line diagnostics. 
We take the statistics from ABOO to be representative of 



(km s-‘) 

Figure 9. Maximum blue-shifted absorption velocities seen in He 
I A10830 vs Ha and H/3 for objects with blue shifted absorption at 
both He I A10830 and either Ha or H/3. The He I A10830 and the 
Balmer lines generally trace the same velocity of outflowing ma¬ 
terial, although there is significant scatter. The H/3 line velocities 
match the He I A10830 velocities less well than the Ha velocities. 
This may be an optical depth effect: the optical depth at H/3 may 
simply be too low to trace the same wind extent as He I A10830 
and Ha, resulting, in general, in smaller observed velocities at H/3. 
Note that the data for the optical lines and the He I A10830 lines 
are not simultaneous. 


CTTSs in general. ABOO examine the following accretion 
and outflow diagnostics: Ha, H/3, Hy, H5, He I A5876, 
Na I D, Ca H K, Ca H A8662, Ca H A8498, Fe I A6I92, 
A5497, and Fe II A4924 and A4352. Our study is well 
suited for comparison with ABOO since we examine all 
of these except for Fe I A6I92. The resolving power of 
their observations {R ^48,000) is also similar to ours in¬ 
creasing the chance that no weak features will be missed 
in one data set but not in the other. We will also com¬ 
pare certain profile statistics from our sample to those of 
Beristain et al. (2001) and Batalha et al. (1996), studies 
which analyzed a smaller number of mass flow diagnos¬ 
tics in similar sample sizes of CTTSs to that of ABOO. 

ABOO find high fractions of both red-shifted absorption 
(12 objects, 40%lg) and blue-shifted absorption (24 ob¬ 
jects, 80%/'/g) in at least one line for the CTTSs in their 
sample. Our sample displays much lower incidences of 
both red and blue-shifted absorption: 26% and 37%, re¬ 
spectively. Contingency tests comparing our study to the 
of Alencar & Basri (2000) are shown in Table 7. The red- 
shifted absorption comparison between our final sample 
and that of ABOO suggests a marginal difference (p=0.I6) 
in the incidence of red-shifted absorption. The blue- 
shifted absorption statistics are signihcantly different for 
both sample comparisons, with p <0.001. 
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Table 7 

Optical contingency table results 


Feature 

Groups 

p-value 

Red absorption 
Blue absorption 

HAEBES vs CTTSs 
HAEBES vs CTTSs 

0.160 

9x10-* 


In addition to the total mass flow statistics, there are 
differences in the absorption statistics in the lines that 
are common to both our study and that of ABOO. This 
comparison is presented in Table 8 and is shown graph¬ 
ically in Figure 10. The most common line displaying 
blue-shifted absorption in their data is Ho followed by 
H/3, Hy, Hh and the Na I D lines. A similar trend is ob¬ 
served in our data except that H/3 shows a slightly higher 
incidence of blue-shifted absorption than Ho. The most 
common red-shifted absorptions features appear at He 
I A5876 and Na I D in their sample. The highest red- 
shifted absorption occurrence in our data is also at He 
I A5876. However, we observe only 1 object with red- 
shifted absorption at Na I D. We find a small incidence 
of inverse P-Cygni profiles in Ha and H/3 (6% and 5%, 
respectively) while ABOO observe a larger incidence of 
red-shifted absorption components at Ha (23%) and H/3 
(23%). ABOO find Fe H A4924 in emission in 40% of their 
sample. Our data show Fe H A4924 emission in only 9% 
of the sample. We also observe both blue and red-shifted 
absorption at Fe H A4352 in a non-negligible fraction 
(13%) of our sample. ABOO do not fit for the absorption 
components at Fe H A4352, likely due to strong nearby 
photospheric features. By comparing the CTTS A4352 
lines in their sample to that of V410 Tau, a WTTS also 


from their sample, we searched for red and blue-shifted 
absorption in the line. We find no profiles that are not 
consistent with purely photospheric absorption. Thus we 
estimate the mass flow incidence in this line to be ~0 in 
the ABOO CTTS sample. 

The largest incidence differences between our sample 
and ABOO, outside of the Balmer lines, occur in the Na 
I D lines. He I A5876, Fe I A5498, and Fe H A4924. At 
He I A5876 we find a higher occurrence of red-shfited ab¬ 
sorption (15%) than ABOO (10%). Our observed fraction 
of He I A5876 red-shifted absorption is more similar to 
that found by Beristain et al. (2001) (16%) for a sample 
of 31 CTTSs. We find significantly lower fractions of red- 
shifted absorption in both Na I D lines (8%) compared to 
ABOO (27%). We find a red-shifted absorption fraction 
of 8% at Fe H A4924. ABOO do not observe any red- 
shifted absorption at Fe H A4924. Similar discrepancies 
are observed in the blue-shifted absorption statistics. At 
Fe I A5498, we find no evidence of blue-shifted absorp¬ 
tion in any objects while ABOO report a detected fraction 
of 60%. Our detected blue-shifted absorption fraction at 
Na I D (^10%) is also much lower than ABOO, ’s (40%). 
Finally, the blue and red-shifted absorption fractions at 
Ha, H/3, Hy, and Hi5 are significantly higher in ABOO’s 
sample of CTTSs than in our HAEBE sample. 

Although we do not focus on the pure emission profiles, 
we note that our detected fraction of emission at He I 
A5876 is low (14%) compared to that of Beristain et al. 
(2001) who detect emission in 100% of their 31 CTTSs. 
In §6 we will comment on possible causes for this differing 
line profile statistic. The possible physical causes of the 
absorption statistics detailed above will be outlined in 
the next section. 


Table 8 

Comparison to line absorption statistics from Alencar & Basri (2000) 




Ha 

(1) 

H/3 

(2) 

H7 

(3) 

H<S 

(4) 

He I 
A5876 
(5) 

Na I 
A5889 
(6) 

Na I 
A5896 
(7) 

Ca H 
A8662 
(8) 

Ca H 
A8498 

(9) 

Ca H 
A3969 
(10) 

Fe I 
A5498 
(11) 

Fe H 
A4924 
(12) 

Fe H 
A4352 
(13) 

ABOO 

t^RA 

9 

8 

8 

8 

3 

8 

8 

4 

1 

1 

0 

0 

0 


% RA 

30 

31 

31 

35 

10 

27 

27 

13 

3 

4 

0 

0 

0 


Nr A 

23 

14 

8 

6 

0 

12 

12 

0 

1 

2 

18 

0 

0 


% BA 

77 

54 

31 

26 

0 

40 

40 

0 

3 

7 

60 

0 

0 

This study 

t^RA 

5 

7 

3 

3 

12 

6 

6 

2 

1 

0 

0 

5 

0 


% RA 

6 

9 

5 

6 

15 

8 

8 

3 

4 

0 

0 

8 

0 


Nr A 

21 

24 

12 

7 

1 

8 

7 

2 

2 

1 

0 

7 

3 


% BA 

24 

32 

20 

14 

1 

10 

9 

3 

8 

2 

0 

11 

5 


6 . DISCUSSION 

The optical line statistics presented in §5 indicate that 
there are differences in the characteristics of, and/or 
the physical mechanisms producing, the mass flows in 
CTTSs vs HAEBES. Below we discuss possible scenarios 
capable of producing the trends and differences seen in 
the measured statistics. 

6 .1. Fe II: Higher ionization states tracing mass flows 
in HAEBES 

Table 8 shows clear evidence that HAEBES, in gen¬ 
eral, show significantly more mass flow activity in the 


Fe H lines than CTTS indicating an absorbing region 
that is hotter or is subject to higher amounts of ionizing 
radition. Fe I A5498 shows a 60% incidence of blue- 
shifted absorption in ABOO’s sample but is not detected 
in absorption in any of our HAEBES. The reverse is true 
for the Fe H A4924 and A4352 lines: we see non-zero 
fractions of both lines exhibiting blue-shifted absorption 
while ABOO do not detect absorption in either of these 
lines in any objects. All of the objects in which we detect 
blue-shifted absorption in Fe H are B-type or early A- 
type objects. This behavior was also noted by Hamann 
& Persson (1992) concerning the relative strengths of Fe 
I and Fe H emission lines in HAEBES versus CTTSs. 
The first ionization potential of Fe is 7.9 eV which corre- 
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Figure 10. Line-by-line mass flow fraction comparison of CTTS and HAEBE samples. The fractions are listed as percentages in Table 8. 
The CTTS data is from ABOO. The CTTSs, in general, show a higher incidence of both blue and red-shifted absorption in most line 
diagnostics compared to HAEBES. The exceptions are at He I A5876 and the Fe II lines. 


spends to a wavelength of 1570 A. The ionizing radiation 
flux for a 10,000 K blackbody, a similar temperature to 
for an AO star, is ~10® times higher at 1570 A 
than for a 3,800 K blackbody, a typical Tgff for a CTTS 
of spectral type MO. In addition, CTTSs display large 
UV excess fluxes due to accretion which can be 10"*^-10^ 
times higher than the stellar photosphere at typical tem¬ 
peratures used to model the excess continuum. Thus the 
actual flux ratio at 1570 A is likely much lower than 10® 
but still favors HAEBES over CTTSs. The larger num¬ 
ber of ionizing photons emitted by B and early A stars 
compared to CTTSs could account for the observed dis¬ 
parities in the Fe II line profile statistics observed here. 

This trend also holds for another low ionization po¬ 
tential metal, Na I, which is ionized by photons with 
energies higher than 5.1 eV, or A<2431 A. The overall 
mass flow incidence in our sample is much lower at Na 
I D than in ABOO’s sample. However, we do observe Na 
I outflow signatures in a majority of our objects that 
also display blue-shifted absorption at Fe H, indicating 
that significant amounts of neutral material are simulta¬ 
neously present in these outflows. Although it is unclear 
to what extent the temperature of the outflow determines 
the ionization state of the wind compared to the ionizing 
flux from the star, Na I is likely not present in the re¬ 
gions of the outflow being sampled by Fe H. The observed 
statistics for Ca H are not significant enough to warrant 
a comparison to those from ABOO. It appears that the 


differences in Fe I and Fe H and Na I morphology statis¬ 
tics between HAEBES and CTTSs suggest more highly 
ionized environments around HAEBES, likely due to the 
higher number of ionizing photons emitted by HAEBES. 

6 .2. Infall velocities 

In Paper I we showed that, on average, the maximum 
red-shifted absorption velocities in HAEBES are smaller 
fractions of the stellar escape velocity than in CTTSs (see 
Figures 5 and 6 from Paper I). The comparison in Paper 
I was performed using only the line velocities measured 
from He I A10830 data. Here we present the same com¬ 
parison in Figure 11 but incorporating all of the optical 
diagnostics displaying red-shifted absorption for each ob¬ 
ject. For HAEBES that show red-shifted absorption in 
more than one diagnostic, their average value is used. We 
note that, as in Paper I, Figure 11 is identical in format 
to Figure 5 from Fischer et al. (2008) in order to facil¬ 
itate easy comparison with the red-shifted He I A10830 
absorption velocities presented in that study. The diago¬ 
nal lines mark the velocity ratios that correspond to the 
infall distance from the stellar surface, given in stellar 
radii, marked above each line. In other words, the dot¬ 
ted line shows the maximum free-fall velocity that can 
be obtained if infall begins at the specified distance from 
the star. We estimate the red-shifted absorption veloc¬ 
ities of 13 CTTSs from the published optical profiles of 
Edwards et al. (1994). Edwards et al. (1994) observed 
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Figure 11. Observed maximum red—shifted absorption velocities 
versus stellar escape velocity for the HAEBES (filled blue circles) 
from this study and CTTSs from Edwards et al. (1994) (filled red 
diamonds). The unfilled plot symbols are the HAEBE (light blue) 
and CTTS (light red) objects from the He I A10830 samples of 
Paper I and Fischer et al. (2008), respectively. Objects common 
to both the optical and He I A10830 samples are connected by a 
dotted line. The dashed lines show, for a given Vesc, the maximum 
velocity obtained if infall begins from the distance labeled next to 
the line (in stellar radii). The HAEBES, on average, show maxi¬ 
mum red—shifted absorption velocities that are smaller fractions of 
their stellar escape velocities, indicating that the infalling material 
originates closer to the stellar surface. 

Balmer line profiles as well as Na I D and He A5876. 
Most of the red-shifted absorption profiles in our sample 
are seen in the same lines so the comparison should not 
be subject to lines tracing mass flows with significantly 
different physical conditions. 

In Figure 11 the CTTSs are plotted as solid red di¬ 
amonds; the HAEBES are represented by the filled 
blue circles. For comparison we have overplotted the 
HAEBES He I A10830 velocities from Paper I as light 
blue unfilled circles and the He I A10830 CTTS velocities 
from Fischer et al. (2008) as unfilled light red diamonds. 
Objects that are common to both the optical and He 
I A10830 samples have their velocities, as measured at 
He I A10830 and in the optical diagnostics, connected 
by a dotted line^. The mean velocity ratio (Ked/14sc) 
for the HAEBE sample is 0.40; for the CTTSs it is 0.78. 
A two-sided KS-test performed on the sample returns 
a p-value of 2.2x10“^ indicating a significant difference 
between the two populations. 

In order to test the effect of uncertainties in determin¬ 
ing Vesc on the KS statistic, we generated ten thousand 
random velocity distributions and performed a KS test 
on each one. We follow the same procedure for assign¬ 
ing uncertainties to the stellar mass and radius (50% and 
25%, respectively) as performed in Paper I. In addition, 
we let the measured red-shifted absorption velocity vary 
by 50%. This provides some accounting of accretion vari¬ 
ability. The result of the KS-tests is shown in the bottom 

^ Two CTTSs are plotted with two different estimates of their 
escape velocities due to different M* and ij, estimates from the 
studies that the spectral line velocities are taken from. 
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Figure 12. Histogram showing the distribution of KS test p- 
values (bin width of 0.0005) for ten thousand trials (bottom panel) 
and EDFs for the HAEBE and CTTS samples (top panel). The 
median value of Pks is marked by the vertical dashed line in the 
bottom panel. In the top panel, the HAEBE EDF is shown in blue; 
the CTTS EDF is shown in red. A majority of the tests (74%) 
return p<0.05 indicating a significant difference between the two 
EDFs. 

panel of Figure 12 where the median p-value is seen to 
be ~0.014. The emipirical distribution function (EDF) 
for each sample is shown in the upper panel. For each of 
the trials each object is assigned a random value within 
the limits of the uncertainty in the object’s escape veloc¬ 
ity and red-shifted absoprtion velocity. The low p-values 
suggest a real difference in the distributions which in turn 
indicates that accretion flows in HAEBES begin deeper 
in the system’s gravitational potential, i.e. closer to the 
star. This supports the similar result found in Paper I, 
although the p-values are, on average, higher in the opti¬ 
cal case due to the red-shifted absorption velocity being 
allowed to vary in the simulation. 

Eigure 11 and Eigure 12 provide additional support 
for our suggestion in Paper I that HAEBES have smaller 
magnetospheres than CTTSs. The evidence in the op¬ 
tical is perhaps even more convincing due to the larger 
number of red-shifted absorption lines used in the com- 
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parison, although there is a large spread in the mea¬ 
sured line absorption velocities for some individual ob¬ 
jects. However, the lack of large p-values in our KS 
test simulation indicates that the velocity ratio differ¬ 
ence between HAEBES and CTTSs is independent of 
the particular line velocity chosen for the comparison. 
The statistical tests also support the idea that the red- 
shifted absorption in HAEBES originates deeper within 
the gravitational potential of the central star. If the rel¬ 
ative distances are similar at which the magnetospheres 
in HAEBES and CTTSs are dominating the kinematics 
of the accretion flow, the velocity ratios of the HAEBE 
and CTTS samples should be comparable. Since they 
are not, there must be some difference in the formation 
mechanisms of the absorption signatures. We attribute 
this difference to smaller magnetospheres in HAEBES 
compared to CTTSs. We note that smaller corotation 
radii, and thus smaller magnetospheres, are suggested for 
most HAEBES compared to CTTSs due to their large 
usini values (see §2 of Paper I; Muzerolle et al. 2004). 
Thus, this result is expected from magnetospheric accre¬ 
tion physics. 

6.3. HAe stars vs HBe stars 

In Paper I we demonstrated the lack of red-shifted ab¬ 
sorption at He I AI0830 in HBe stars (Figure 7 of Paper 
I). This same trend holds for our larger optical sample. 
This is shown in Figure 13, which plots the difference 
between the blue (v < 0) and red {v > 0) line equivalent 
widths versus Te//, for the objects with at least one pro¬ 
file displaying evidence of blue or red-shifted absorption. 
If an object shows more than one profile in absorption, we 
average the equivalent width difference. The equivalent 
width difference is essentially a measure of profile type: 
red-shifted absorption results in a negative equivalent 
width difference while blue-shifted absorption results in 
a positive equivalent width difference. Objects in each 
quadrant are color-coded for clarity. There is an absence 
of objects in the upper-left quadrant. Although V1818 
Ori shows evidence of red-shifted absorption in the Na I 
D lines, we do not calculate an equivalent width for ei¬ 
ther line due to the strong interstellar absorption features 
superimposed on the line. Thus there are only 3 HBe ob¬ 
jects with red-shifted absorption shown in Figure 13 and 
they all lie very close to the vertical dashed line. The IPC 
and red absorption He I A5876 profiles generally have 
small equivalent width differences (0.0 A>H4-Wi.>-2.0 
A) and comprise most of the red-shifted absorption ob¬ 
jects (red crosses and magenta right-facing triangles). 

This result is also born out in the overall profile statis¬ 
tics: 15 of the 19 red-shifted absorption objects are HAe 
stars; 4 are HBe stars, with 2 of those having spectral 
types of B9. There are 45 HAe stars in the final sam¬ 
ple and 33 HBe stars. A contingency test comparing 
HAe and HBe red-shifted absorption returns p=0.04 pro¬ 
viding support for a real difference in the morphology 
statistics. A test comparing the blue-shifted absorp¬ 
tion reveals no difference between HAe and HBe objects 
(p=0.48). A similar trend is seen in the He I A10830 data 
but the larger optical sample has added significance to 
the red-shifted absorption comparison. 

The similarity in the optical and He I A10830 HAe 
versus HBe comparisons suggests that the conclusions 



Figure 13. The T^ff vs blue equivalent width minus red equiva¬ 
lent width for all of the objects in our sample that have absorption 
profiles. Objects with more than one line in a quadrant are shown 
using the average of the equivalent width differences for those lines. 
Colors are used for clarity in distinguishing quadrants. The trend 
found for the He I A10830 data from Paper I holds for the optical 
sample: there is a lack of red—shifted absorption profiles in HBe 
stars. The strongest blue-shifted absorption profiles are found in 
the HBe group. 

reached in Paper I are not altered by the analysis pre- 
sented here. HAe stars exhibit red-shifted absorption 
more frequently than HBe stars indicating that mass in¬ 
fall, most likely magnetically channeled, is more common 
in HAe stars. Outflows seem to be present in similar frac¬ 
tions of HAe and HBe objects suggesting that the wind 
mechanism is equally as efficient in both groups of ob¬ 
jects. 

6.4. Simultaneous blue and red-shifted absorption and 
accretion driven outflows 

In Paper I we noted the lack of simultaneous blue and 
red-shifted absorption in the He I A10830 profiles. To 
clarify, simultaneous absorption signatures in the cur¬ 
rent, optical context refer to detected absorption in any 
diagnostics from the same exposure, or exposures taken 
very close in time to one another. This includes simulta¬ 
neous red and blue-shifted detections in the same line, 
such as that seen in He I A10830 for the CTTSs from Ed¬ 
wards et al. (2006) , or detections in two different lines as 
seen in Alencar & Basri (2000). This lack of simultane¬ 
ous absorption contrasted sharply with the large detected 
fraction of simultaneous absorption signatures in CTTSs 
presented by Edwards et al. (2006). 

The optical data continues this trend: we observe only 
1 object (HD 163296) that shows clear red and blue- 
shifted absorptions signatures in at least one line from 
the same exposure. Based on the published centroid ve¬ 
locities of the detected absorption features from ABOO, 
we count 17 out of 30 CTTSs, or 57%, showing red and 
blue-shifted absorption in at least one line from nearly 
simultaneous exposures. The lack of simultaneous ab¬ 
sorption features observed in our optical sample supports 
the arguments presented in Paper I, namely that accre- 
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tion driven outflows in HAEBES may be more difficult 
to produce. The strong optical outflow signatures ob¬ 
served in both HAe and HBe objects that do not also 
exhibit red-shifted absorption, combined with the con¬ 
clusion from Paper I that the He I A10830 PC morpholo¬ 
gies are likely accretion-driven, suggests that accretion 
in some HAEBES may proceed through a non-magnetic 
mechanism such as a boundary layer. This would explain 
the relative lack of simultaneous red and blue-shifted ab¬ 
sorption signatures in HAEBES compared to CTTSs. 

It is still unclear whether or not the accretion flows 
detected in the red-shifted absorption objects can drive 
outflows. Outflow signatures in these objects are not 
observed (except for HD 163296) suggesting that if the 
accretion is related to driving outflows, the outflow sig¬ 
natures are too weak to be detected or are launched in 
a geometry that requires a specific viewing angle to ob¬ 
serve (e.g. a highly collimated jet emerging from the 
polar region). This seems unlikely since these simultane¬ 
ous signatures are almost entirely absent from our sam¬ 
ple, suggesting that they are present in very few objects. 
Indeed, jets have been detected in a relatively fewer num¬ 
ber of HAEBES (e.g. Grady et al. 2004; Ellerbroek et al. 
2014) compared to CTTSs hinting that they do not form 
as efficiently in HAEBES. 

One potential explanation is that, in general, the ac¬ 
creting material as traced by red-shifted absorption is 
not energetic enough to provide the energy necessary 
to drive an outflow. Support for this idea is seen in 
Figure 11 from this study and Figure 5 from Paper I 
where the red-shifted absorption velocities of HAEBES 
are found to be smaller fractions of their stellar escape ve¬ 
locities. We also find support for this scenario in Figure 
7 from Paper I which demonstrates the relative lack of 
objects with low accretion rates that also exhibit blue- 
shifted absorption profiles. We note that this distinc¬ 
tion is not clear in our optical data, i.e. the separation 
of profile type with accretion rate is not distinct. We 
have fewer estimates of M in hand, however, for the red- 
shifted absorption objects from the optical sample and 
thus the comparison is less significant. This could be im¬ 
proved and further tested by obtaining simultaneous line 
profile morphologies and accretion rate estimates for a 
large sample of HAEBES. 

The idea of accretion powered stellar winds was first 
proposed for CTTSs by Matt & Pudritz (2005). They 
suggest that a portion of the deposited accretion energy 
will heat the corona (via conduction or magnetosonic 
wave dissipation) and assist in launching stellar winds. 
While it is unclear whether or not this process is oc¬ 
curring in CTTSs, magnetically channeled accretion, in¬ 
volving smaller magnetospheres in HAEBES compared 
to CTTSs, will provide a smaller fraction of the nec¬ 
essary energy for launching an outflow from the stellar 
surface of HAEBES than from the surface of CTTSs. 
Since the outflows observed in CTTSs are generally ac¬ 
cepted to be driven largely by accretion (e.g. Hartigan 
et al. 1995; Edwards et al. 2006), the amount of energy 
deposited by the accretion flow onto the stellar surface 
could strongly influence the strength of the outflow. For 
HAEBES, if we imagine that all of the kinetic energy 
of the accretion flow is converted into kinetic energy in 
a wind driven from very near the stellar surface (i.e. a 


stellar wind), these outflows will be incapable of escaping 
the gravitational potential of the star due to the small 
(^0.5; Figure 11) fractions of the stellar escape velocity 
that the accretion flows attain during their infall. Ob¬ 
viously this is an ideal scenario, but the energy budget 
argument holds if a large portion of the energy in the 
wind is provided by the accretion flow, which it seems 
to be in CTTSs. To summarize, the lower infall veloci¬ 
ties relative to the stellar escape velocities in HAEBES 
may indicate that these accretion flows cannot efficiently 
drive mass outflows from the stellar surface. 

If magnetically controlled mass infall in HAEBES is 
incapable of driving outflows, another form of accretion 
may be partly responsible for driving the outflows from 
objects with large M. We suggested in Paper I that 
boundary layer accretion in HAEBES could provide the 
necessary energy, although it is unclear exactly how the 
energy is tranferred to the outflow. One possible sce¬ 
nario discussed in Paper I for launching outflows from 
the boundary layer is the Shu et al. (1988) X-celerator 
mechanism. The high vs'mi values of HAEBES are ideal 
for producing this type of outflow. Better M estimates 
for a larger number of HAEBES will help elucidate the 
relationship between profile type and the role of the ac¬ 
creting material in producing the outflow signatures. 

6.5. He I X5876 emission lines 

Beristain et al. (2001) examine the He I A5876 in a 
sample of 31 CTTSs. They find an emission component 
(broad, narrow, or both) in all of the objects in their 
sample. They interpret the broad component emission 
as arising in a hot wind as well as in the accretion flow 
onto the star. The narrow emission component is inter¬ 
preted as forming in the post-shock gas at the base of 
an accretion column at the surface of the star. Our sam¬ 
ple contains 13 HAEBES that display emission profiles 
at He I A5876. None of these objects show the narrow 
component that is commonly seen in CTTSs, as well as 
WTTSs. This may be due to the much higher average 
nsini values of HAEBES compared to CTTSs. However, 
11 of the 13 objects show broad emission profiles similar 
in width to those observed by Beristain et al. (2001) and 
8 objects show unambiguous blue-shifted absorption in 
a separate but simultaneous diagnostic. Of these 8 ob¬ 
jects, five show unambiguous evidence of mass outflow 
at He I A10830, although these observations are not si¬ 
multaneous to the He I A5876 observations. Due to the 
prevalence of blue-shifted absorption observed in the He 
I A5876 emission objects, our data supports the sugges¬ 
tion by Beristain et al. (2001) that the broad He I A5876 
emission component may partly arise in a hot wind near 
the star. 

7. SUMMARY AND CONCLUSIONS 

We have presented a large sample census of mass flows 
around Herbig Ae/Be stars using a wide variety of mass 
flow diagnostics and determined the fraction of objects 
that display unambiguous (i.e. red and blue-shifted ab¬ 
sorption) mass flow signatures. We find that the Na I 
D interstellar absorption line centers are good approxi¬ 
mations, to within 8 . 1 j'j 3 g km s“^, to the center-of-mass 
radial velocity in most HAEBE systems, confirming the 
conclusion reached by Finkenzeller & Jankovics (1984). 


24 


Of the 78 HAEBES included in the final analysis, 26%^5 
show direct signs of mass infall and 37%1| show direct 
evidence of mass outflow. These percentages are very 
similar to those found in Paper I using the He I A10830 
line. 

Overall, the optical data reinforce the conclusions out¬ 
lined in Paper I. More specifically, we find lower rates of 
both red and blue-shifted absorption for HAEBES than 
is observed for CTTSs, suggesting that, in general, the in¬ 
nermost environments of HAEBES are not simply scaled 
analogues of CTTSs. Similar to the He I A10830 data, 
the HAEBE optical lines show maximum red-shifted ab¬ 
sorption velocities that are smaller fractions of the stel¬ 
lar escape velocities than is found for CTTSs. This adds 
support to the argument presented in Paper I concerning 
the smaller size of HAEBE magnetospheres compared to 
those of CTTSs. This has implications for the ability 
of HAEBE accretion flows to power strong winds: less 
energetic funnel flows provide less of the required energy 
for launching outflows. This is supported by the lack of 
simultaneous blue and red-shifted absorption signatures 
seen in both the optical sample and the He I A10830 sam¬ 
ple from Paper I. This contrasts with CTTSs for which 
the accretion process is believed to be the main energy 
source for outflows. 

We also find a significant difference between the occur¬ 
rence of red-shifted absorption in HAe stars compared to 
HBe stars; both groups appear to exhibit outflow signa¬ 
tures at a similar rate. This suggests that HBe stars ac¬ 
crete material from their disks differently than HAe stars, 
perhaps through a boundary layer instead of along mag¬ 
netic field lines. The prominence of blue-shifted absorp¬ 
tion in accreting objects in both subgroups suggests that 
the outflows in HAEBES are accretion driven. However, 
the lack of simultaneous red and blue-shifted absorption 
in all but one HAEBE (HD 163296) indicates that what¬ 
ever geometry is producing the absorption signature is 
not capable of producing the opposite absorption sig¬ 
nature. In other words, boundary layer accretion, which 
should not produce red-shifted absorption, may be capa¬ 
ble of driving strong outflows which results in observed 
blue-shifted absorption; magnetically controlled accre¬ 
tion from a small magnetosphere produces red-shifted 
absorption but is perhaps not energetic enough to drive 
outflows from the stellar surface, resulting in the lack 
of simultaneous red-shifted absorption in these profiles. 
Due to the observed statistics, this suggests a general 
transition from magnetically controlled accretion in HAe 
stars to boundary layer accretion in HBe stars. 
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APPENDIX 


OPTICAL LINE PROFILES 

Below are the extracted profiles for a piece of the sample. The profiles for the entire sample will be available with 
the published version. The format of each figure is the same as Figure 3 and Figure 4. 
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Figure 14. Extracted line profiles for AB Aur. 
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Figure 15. Extracted line profiles for AE Lep. 
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Figure 16. Extracted line profiles for BD+61 154. 
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Figure 17. Extracted line profiles for BF Ori. 
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Figure 18. Extracted line profiles for BH Cep. 
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Figure 20. Extracted line profiles for DW CMa. 
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Figure 23. Extracted line profiles for HD 142666. 
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Figure 24. Extracted line profiles for HD 163296. 
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Figure 26. Extracted line profiles for HD 190073. 










































































MgI5167 


MgI5173 


MgI5184 


iBbBII 


NO DATA 


[Nil] 6548 


[Nil] 6583 


[Nil] 5755 




[OI] 5577 


[01] 6300 


[OI] 6363 


[OIII] 5007 


■MBMMliMM 


-400 -200 0 200 400 -400 -200 0 200 400 -400 -200 0 200 400 PC IPC DP 


Velocity (km s"') 


Figure 27. Extracted line profiles for HD 203024. 
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Figure 28. Extracted line profiles for HD 244314. 
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Figure 29. Extracted line profiles for HD 244604. 
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Figure 31. Extracted line profiles for HD 249879. 





























































